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ABSTRACT 


A  moving-coordinate  prediction-model  technique  applied  to  East  Coast 
cyclones  has  been  extended  to  predict  cyclone  behavior  in  most  of  North  America 
for  both  summer  and  winter.  Equations  based  on  stratification,  incorporation  of 
past  history,  and  addition  of  derived  variables  are  developed  from  833  cases. 
Comparison  with  equations  based  on  unstratified  data,  no  past  history,  and  only 
simple  variables  shows  that  the  latter  set  of  equations  is  usually  superior.  This 
“base-technique”  set  of  equations  yields  vector-position  rms  errors  of  2.38,  4.26, 
and  5.96  deg  lat  and  central-pressure  rms  errors  of  4.28,  6.88,  and  9.03  mb  for 
12,  24,  and  36  hr  when  applied  to  213  independent  winter  cases.  Subsequent  refine¬ 
ments  achieve  vector-position  rms  errors  of  2.31,  4.10,  and  5.88  deg  lat  and 
central-pressure  rms  errors  of  4.19,  6.88,  and  9.03  mb  for  the  3  time  periods. 
Comparable  results  for  climatology  are  3.42,  6.13,  and  8.34  deg  lat  for  position 
and  4.83,  8.03,  and  10.50  mb  for  central  pressure.  Results  for  summer  cyclones 
are  also  presented. 

Equations  based  on  surface  data  only  are  developed  and  tested.  They  are 
useful  whenever  upper-air  data  are  unavailable. 

Worksheets  for  operational  application  of  the  technique  are  included. 
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1.0  INTRODUCTION 


This  study  represents  an  extension  of  previous  work  dealing  with  the  prediction 
of  cyclones  over  the  eastern  United  States  [14]  and  cyclones  and  anticyclones  over 
Europe  and  Asia  [7] . 

The  method  for  predicting  the  movement  and  intensification  of  a  cyclone 
consists  of  prediction  equations  derived  statistically  and  based  on  a  moving- 
coordinate  system  (a  system  in  which  predictor  information  is  measured  at  points 
fixed  with  respect  to  the  moving  cyclone  rather  than  at  points  fixed  with  respect  to 
the  earth). 

Successful  results  on  East  Coast  cyclones  [14]  (which  showed  the  method  to 
be  competitive  with  the  skill  of  synoptic  forecasters),  coupled  with  the  relative 
simplicity  with  which  a  forecast  could  be  completed  on  just  a  desk  calculator, 
made  the  method  feasible  for  operational  application.  During  the  past  few  winters, 
the  prediction  equations  for  East  Coast  cyclones  have  been  used  either  operationally 
or  experimentally  as  a  forecast  aid  at  several  East  Coast  weather  stations,  including 
the  National  Weather  Analysis  Center  (NWAC).  This  apparent  acceptance  of  the 
method  warranted  extension  of  the  study  to  derive  prediction  equations  for  other 
geographical  areas. 

The  prediction  of  cyclone  displacement  and  intensification  is  not  an  end  in 
itself:  it  should  provide  information  helpful  in  preparing  a  forecast  of  the  associated 
weather. 

The  system  of  equations  evolving  from  this  study  can  be  used  operationally 
either  manually  with  the  aid  of  a  desk  calculator  or  completely  automatically  with 
the  aid  of  a  computer. 
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2.0  DEFINITION  OF  THE  PROBLEM 


The  objective  of  this  study  was  to  derive  equations  for  predicting  the  movement 
and  intensification  of  surface  cyclones  over  an  area  of  North  America  for  intervals 
of  from  12  to  36  hr  in  advance. 

2.1  Area  Studied 

The  area  defined  in  Fig.  2-1,  which  includes  most  of  North  America,  was 
chosen  for  this  study.  Previous  similar  studies  involved  East  Coast  (United  States) 
cyclones  [14J  and  European  and  Asian  cyclones  and  anticyclones  [7] . 

2.2  Selection  of  Cases 

Table  2-1  shows  the  number  of  cyclones  selected  for  each  of  the  zones  of 
Fig.  2-1.  The  cyclones  were  selected  by  examining  all  0000-  and  1200-GCT  surface 
charts*  for  the  winters  (November— March)  of  1955-1956  through  1959-1960  and 
summers  (May— September)  of  1955  through  1959.  A  cyclone  was  accepted  if  it 
retained  its  identity  for  at  least  36  hr  in  winter  or  24  hr  in  summer. 

The  sample  of  cyclones  selected  for  the  winters  of  1955-1956  through  1958- 
1959  was  designated  the  dependent  sample  for  derivation  of  winter-cyclone  equations; 
the  sample  of  summer  cyclones  for  the  summers  of  1955  through  1958  was  designated 
the  dependent  sample  for  derivation  of  summer-cyclone  equations.  Independent 
samples  were  made  up  of  1959-1960  winter  cyclones  and  1959  summer  cyclones. 

♦Microfilm  copies  of  surface  maps,  analyzed  by  the  National  Weather  Analysis 
Center,  were  used.  From  these,  it  was  possible  to  construct  maps  showing  individ¬ 
ual  cyclone  tracks  at  12-hr  intervals,  which  served  as  the  basis  for  selection. 
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Fig.  2-1.  The  area  for  which  cyclones  were  studied.  The  ai’ea  is  divided  into  four  zones  (northwest,  north¬ 
east,  southwest,  and  southeast)  by  latitude  458N  and  longitude  90°\V. 


TABLE  2-1 

NORTH  AMERICAN  CYCLONES  SELECTED  FOR  STUDY, 
1955-1960 


Zone* 

No.  of 

winter  eye lones 

No.  of  summer  cyclones 

1955-1959t 

1959-1960$ 

Total 

1955-1958t 

HB^SI 

Tota  1 

NW 

1> 

43 

217 

226 

79 

305 

NE 

15 4 

34 

188 

199 

38 

237 

SW 

220 

69 

289 

129 

20 

1 49 

SE 

285 

67 

352 

1 46 

32 

178 

A  1  1 

833 

213 

io46 

700 

169 

869 

♦See  Fig.  2-1. 
f Dependent  sample, 
j Independent  sample. 
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3.0  THE  PREDICTION  TECHNIQUE 
3.1  Screening  Regression 

Multiple-linear-regression  analysis  was  used  to  relate  cyclone  movement 
and  intensification  (predictands)  with  certain  atmospheric  variables  (predictors) 
thought  to  have  a  bearing  on  cyclone  behavior.  The  literature  shows  that  many 
predictors  have  been  suggested.  Although  it  would  seem  desirable  to  utilize  all 
predictors  deemed  important  on  the  basis  of  previous  theoretical,  Bynoptic,  and 
empirical  work,  Lorenz  [5]  points  out  that  a  prediction  equation  should  contain 
few  predictors  in  comparison  with  the  size  of  the  developmental  sample;  if  there 
are  too  many,  a  relationship  that  fits  the  sample  used  to  establish  it  is  likely  to 
fail  when  applied  to  a  new  sample.  The  screening  procedure  of  Miller  [6]  was 
used  in  this  study  to  limit  the  number  of  predictors. 

The  object  of  the  screening  procedure  is  to  select  from  a  set  of  possible 
predictors  the  subset  that  contributes  most  significantly  and  independently  to 
reducing  the  variance  of  the  predictand.  From  an  array  of  possible  predictors, 
the  screening  procedure  selects  first  the  one  that  has  the  highest  linear  corre¬ 
lation  with  the  predictand  in  question.  This  predictor  is  then  held  constant,  and 
partial-correlation  coefficients  between  the  predictand  and  each  of  the  remaining 
predictors  are  examined;  the  predictor  now  associated  with  the  highest  coefficient 
is  the  second  one  selected.  Additional  predictors  are  chosen  similarly  until  a 
selected  predictor  fails  to  explain  a  significant  additional  percentage  of  the 
remaining  variance  of  the  predictand. 

The  criterion  of  significance  as  applied  to  the  screening  procedure  is  not 
clear  cut  because  the  usual  F-test  methods  [e.g.,  12]  are  not  applicable  [6] .  If 
a  predictor  is  chosen  at  random  from  a  group  of  predictors,  an  F-test  is  usually 
taken  at  the  95%  level;  this  allows  for  a  l-in-20  chance  of  considering  the  pre¬ 
dictor  significant  when  in  fact  it  is  not.  Because  the  screening  procedure  does 
not  select  predictors  randomly,  a  more  severe  test  is  needed  to  specify  a  l-in-20 

i 

chance.  For  his  screening  procedure,  Miller  suggested  [6]  that  the  critical  F-value 
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be  a  function  of  the  number  of  possible  predictors.  The  F-test  was  used  in  this 
form  in  this  study.  After  the  significant  predictors  have  been  selected,  the 
regression  coefficients  are  obtained  by  the  method  of  least  squares. 

3.2  The  Moving-coordinate  Grid 

The  grid  for  extracting  predictor  information  accompanies  the  cyclone  as 
it  moves,  so  variables  are  measured  at  locations  fixed  with  respect  to  the  cyclone 
center  rather  than  with  respect  to  specific  geographical  locations.  The  grid  is 
shown  in  Fig.  3-1.  The  gridpoint  defined  by  the  (k,i  )-location  (10,5)  is  placed  at 
the  center  of  the  cyclone,  and  the  grid  is  oriented  so  that  the  line  k  =  10  coincides 
with  the  meridian  passing  through  the  center  of  the  cyclone.  In  practice,  grid 
placement  and  data  tabulation  were  done  by  computer  programs,  and  'Analyzed 
maps"  were  on  magnetic  tape.  On  a  polar  stereographic  projection  with  standard 
parallel  at  60°N,  the  17  x  13  array  forms  a  set  of  evenly  spaced  points.  On  a  map 
scale  of  1:15,000,000,  the  grid  interval  is  1  in.  This  is  the  same  interval  as  in 
the  JNWP  grid  [2] .  At  60°N,  where  the  scale  is  true,  one  grid  interval  equals 
381  km.  The  221  points  defined  by  this  grid  system  were  the  ones  used  for  basic 
predictor  tabulation  by  superimposing  the  grid  on  analyzed  maps  of  sea-level 
pressure,  700-mb  height,  and  500- mb  height. 

3.3  Predictands 

Latitudinal  displacement,  longitudinal  displacement,  and  change  in  central 
pressure  were  chosen  (Table  3-1)  as  12-,  24-,  and  36-hr  predictands  for  winter- 
cyclone  equations  and  as  12-  and  24-hr  predictands  for  summer-cyclone  equations. 
The  same  NWAC  surface  maps  used  to  select  the  cases  (see  Section  2.2)  were  used 
as  the  source  of  predictand  data. 

3.4  Predictors  Considered 

Hemispheric-data  tapes  [11]  were  used  as  the  primary  source  of  predictor 
data.  Special  preprocessing  programs  (see  Appendix  A)  automatically  derived 
gridpoint  arrays  of  pressure,  height,  and  thickness  data  for  each  cyclone  in  the 
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k  - fr. 

Fig.  3-1.  Moving-coordinate  grid  overlay.  Gridpoint  (10,5) 
in  the  (k,f)-array  is  always  positioned  at  the  center  of  the  cyclone, 
and  the  grid  is  always  oriented  so  that  the  line  k  10  coincides 
with  the  meridian  passing  through  the  center  of  the  cyclone.  One 
grid  interval  equals  one  JNWP  grid  interval  (381  km  at  60°N). 
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TABLE  3-1 
PREDICTANDS 


Descr iption 

Symbo 1 

Unit  of 
measurement 

12-hr  northward  displacement 

A 

N 

12 

deg  lat 

12-hr  eastward  displacement 

*12 

deg  lat 

12-hr  change  in  central  pressure 

°12 

mb 

24-hr  northward  displacement 

N  , 

24 

deg  lat 

24-hr  eastward  displacement 

A 

E24 

deg  lat 

24-hr  change  in  central  pressure 

A 

0  . 

24 

mb 

36-hr  northward  displacement* 

R  , 

36 

deg  lat 

36-hr  eastward  displacement* 

deg  lat 

36-hr  change  in  central  pressure* 

*# 

mb 

*For  winter  cyclones  only. 
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TABIC  w 

FJ0S5I ni  E  PREDICTORS 


Present  latitude  of  cyclone 
Present  longitude  of  cyclono 
5oh-level  pressuro 
12-hr  sea- level  pressure  change 


/uo-mb  height 

12-hr  /OO-mb  height  change 
500-mb  height 

12-hr  500-mb  height  change 

1 000- 700-mb  thickness 

12-hr  1000- 700-mb  thickness  change 

1000-500-mb  thickness 

12-hr  1000-500-mb  thickness  change 


500-mb  absolute  vorticity 


12-hr  vorticity  change 

Then*# I  vorticity  (1000-500-mb) 


12-hr  therma I  vorticity  change 


500-mb  vorticity  advection 


1000- 700-mb  thickness  advection 


1C00  500-mb  thickness  advection 


Thermal  vorticity  advection 


lap  loci  an  of  the  1000-500-mb 
thickness  advection 

700-mb  geostrophic  wind 

12-hr  700-mb  geostrophic  wind  change 

500-mb  geostrophic  wind 

12-hr  500-mb  geostrophic  wind  change 

1000-500-mb  thermal  wind 

12-hr  1000-500-mb  thermal  wind  change 

1000- 700-mb  thermal  wind 

12-hr  1000- 700-mb  thermal  wind  change 


Previous  6-hr  latitude  change 
Previous  6-hr  longitude  change 
Previous  6-hr  central-pressure  change 
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developmental  sample.  These  basic  predictors  are  represented  by  types  A  and  B 
in  Table  3-2. 

More  'fcomplex"  predictors  (type  C,  Table  3-2),  formed  from  various  com¬ 
binations  of  the  basic  predictors,  were  generated  on  the  IBM  7090  by  conventional 
finite -difference  procedures.  The  Laplacian  of  the  1000— 500-mb  thickness 
advection  was  suggested  by  Petterssen  [9]  as  important  to  cyclone  development. 
Advection  of  thermal  vorticity  by  the  thermal  wind  is  also  known  as  the  'develop¬ 
ment"  or  "baroclinic"  term  [13].  East-west  and  north-south  components  and  mag¬ 
nitudes  of  the  geostrophic  wind  over  the  cyclone  were  derived  at  700  and  500  mb  to 
be  used  as  measures  of  steering.  Various  steering  techniques  [1,  4,  8,  10]  for 
forecasting  the  displacement  of  cyclones  are  widely  used  by  practicing  forecasters. 
Twelve-hour  changes  in  thermal  wind  components  were  also  computed  to  include 
steering  tendencies  [e.g.,  3]  as  possible  predictors. 

Miller's  screening- regression  technique  (as  programmed  for  the  IBM  7090) 
could  simultaneously  examine  only  175  predictors,  so  the  number  of  possible 
predictors  was  reduced  subjectively  before  screening  regression  was  applied. 
Section  5.0  discusses  the  reduction  in  more  detail. 

3.5  Gross-error  Checking 

An  important  preliminary  step  in  the  preparation  of  the  data  for  the  prediction 
experiments  was  a  systematic  man-machine  search  for  and  correction  of  gross 
errors.  Suspicious  gridpoint  values  of  parameters  were  identified  by  a  series  of 
tests  performed  by  the  IBM  7090,  and  listings  of  these  values  were  compared  with 
microfilm  copies  of  manuscript  maps;  where  necessary,  corrections  were  made  to 
the  data  filed  on  magnetic  tapes.  The  details  of  the  computer  programs,  data- 
handling  procedures,  and  error-checking  procedures  are  given  in  Appendix  A. 
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4.0  SYNOPTIC  CLIMATOLOGY 


Preparatory  to  deriving  prediction  equations  for  cyclones,  their  climatological 
characteristics  and  their  supposedly  related  atmospheric  variables  were  computed 
to  aid  in  suggesting  possible  predictors.  To  determine  the  difference  in  character¬ 
istics  for  different  regions  of  the  area  (Fig.  2-1),  the  means  and  standard  deviations 
of  all  predictands  and  possible  predictors  were  computed  for  the  entire  area 
(unstratified)  and  for  regions  within  the  area  (stratified). 

4.1  Winter  Cyclones 

Table  4-1  contains  the  means  and  standard  deviations  of  the  northward  and 
eastward  displacements  and  changes  in  central  pressure  of  833  North  American 
winter  cyclones.  Figure  4-1  depicts  the  mean  cyclone  tracks  in  each  of  the  four 
zones.  The  southeastern- zone  cyclones  undergo  the  maximum  deepening  (12  mb 
in  36  hr);  the  northwestern  zone  shows  no  change  in  central  pressure  for  36  hr. 
These  two  zones  also  contained  the  faster-moving  cyclones  (about  25  knots);  the 
southwestern-  and  northeastern-zone  cyclones  averaged  about  20  knots.  The 
southeastern  zone,  of  course,  contains  cyclones  developing  along  the  East  Coast, 
where  the  thermal  contrast  between  land  and  water  is  conducive  to  cyclone 
intensification.  The  northwestern  zone  contains  the  cyclones  that  develop  in  south¬ 
western  Canada.  The  characteristic  southeastward  track  of  these  "Alberta" 
cyclones  shows  up  in  the  mean  track  for  northwestern-zone  cyclones. 

The  associated  mean  maps  of  pressure,  height,  and  thickness  are  shown  as 
Figs.  4-2  through  4-7. 

The  differences  among  the  four  sea-level  pressure  maps  (Fig.  4-2)  are  not 
very  spectacular.  The  mean  locations  and  central  pressures  for  the  cyclones  are 
given  in  Table  4-2.  The  northeastern-zone  cyclones  are  generally  the  deepest 
because  this  zone  contains  many  of  the  more  mature  cyclones,  which  have  already 
undergone  most  of  their  development.  Most  of  the  zones  have  an  area  of  high 
pressure  to  the  west  and  northwest  of  the  cyclone.  The  main  exception  is  the 
northwestern  zone,  where  the  high  pressure  to  the  west  is  very  weak  and  another 
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Fig.  4-1.  Mean  tracks  of  winter  cyclones  in  North  America  by  zone,  1955-1959 
(dependent  sample).  0  12-hr  displacement;  A  24-hr  displacement;  •  36-hr  dis¬ 
placement.  Value  adjacent  to  symbol  refers  to  mean  change  in  central  pressure 
(millibars). 
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Fig.  4-3.  Mean  maps  of  12-hr  sea-level  pressure  change  for  winter  cyclones  in  North  America,  1955-1959 
(dependent  sample).  Isallobars  are  labeled  in  millibars. 
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(dependent  sample).  Isallopachs  are  labeled  in  tens  of  feet. 


cyclone  is  evident  in  the  upper-left  corner.  This  is  undoubtedly  due  to  cyclonic 
activity  in  the  Gulf  of  Alaska.  Similarly,  the  low  pressure  in  the  upper-right 
corner  of  the  northeastern  zone  is  due  to  the  Icelandic  Low. 

The  12-hr  pressure -change  maps  (Fig.  4-3)  are  all  similar  in  that  each 
has  a  katallobaric  center  east  of  the  cyclone  and  an  anallobaric  center  west.  The 
magnitude  of  these  centers  is  larger  in  the  two  northern  zones  than  in  the  southern. 

At  500  mb  (Fig.  4-4),  the  cyclone  is  located  under  a  southwesterly  flow  in 
all  zones  except  the  northwestern,  where  the  flow  is  westerly.  This  is  probably 
a  result  of  many  cases  of  cyclones  imbedded  in  a  northwesterly  flow  pattern, 
which  is  not  uncommon  for  that  zone.  In  the  other  zones,  where  the  flow  is  south¬ 
westerly,  the  trough  line  is  located  200  or  300  naut  mi  upstream  from  the  cyclone. 
The  most  pronounced  trough  occurs  in  the  northeastern  zone. 

The  maximum  12-hr  height  falls  at  500  mb  (Fig.  4-5)  are  located  about  200 
naut  mi  west  of  the  cyclone  in  all  four  zones.  The  general  pattern  shows  rise 
centers  east  and  west  of  the  falls  for  each  zone.  The  greatest  fall  center,  in 
excess  of  240  ft,  is  in  the  northeastern  zone.  The  riBe  center  ahead  of  the  cyclone 
is  toward  the  northeast  or  east— except  in  the  northwestern  zone,  where  it  is 
toward  the  southeast.  The  rather  small  magnitudes  of  both  maxima  and  minima 
can  be  accounted  for  only  by  the  apparent  variability  in  location  of  the  rise  and 
fall  centers  among  individual  cases,  tending  to  damp  them  in  the  mean. 

The  orientation  of  the  1000— 500-mb  thickness  contours  (Fig.  4-6)  is 
similar  to  the  500-mb  height  field  (Fig.  4-4).  The  absolute  values  of  the  thick¬ 
ness  contours  are  about  600  ft  lower  in  the  northern  zones  than  the  southern, 
which  is  simply,  a  reflection  of  the  temperature  variation  with  latitude  in  the 
mid-troposphere . 

The  12-hr  thickness  change  for  the  1000— 500-mb  layer  (Fig.  4-7)  is 
proportional  to  the  change  in  mean  temperature  for  the  layer.  It  is  likely  that 
this  temperature  change  provides  a  rough  estimate  of  the  thermal  advection 
for  the  layer.  In  the  mean-thickness-change  maps,  negative  values  refer  to 


20 


sample).  Isobars  are  labeled  in  millibars. 


I 

I 

I 

1 

I 

I 

1 

I 

I 

1 

I 

I 


•»*> 


i 


sample).  Isohypses  are  labeled  in  tens  of  feet. 


I 


I 

1 

I 

1 

l 

T 


i 


26 


cold  advection  and  positive  values  to  warm.  The  relatively  intense  cold  and  warm 
advections  in  the  northeastern  zone  result  from  the  strong  temperature  contrasts 
associated  with  the  well -developed  mature  cyclones  that  typify  the  northeastern 
zone. 


Figures  4-8  through  4-13  are  simply  the  results  of  pooling  the  four  zones 
shown  in  Figs.  4-2  through  4-7  and  are  representative  of  winter  cyclones  in  the 
entire  North  American  cyclone  prediction  area. 

4.2 Sum  me  r  Cyclones 

Table  4-3  contains  the  means  and  standard  deviations  of  the  northward  and 
eastward  displacements  and  changes  in  central  pressure  of  700  North  American 
summer  cyclones  (Table  2-1).  No  stratification  by  geographical  zone  was  per¬ 
formed.  Comparison  of  Table  4-2  with  the  "All11  portion  of  Table  4-1  indicates 
that  summer  cyclones  move  more  slowly  and  deepen  less  than  winter  cyclones. 
Also,  the  standard  deviations  show  that  summer  cyclones  vary  less  in  movement 
and  intensification  than  winter  cyclones. 

Figures  4-14  through  4-19  contain  the  mean  pressure,  height,  and  thickness 
maps  for  the  summer-cyclone  sample.  A  comparison  with  the  mean  maps  of 
winter  cyclones  (Figs.  4-8  through  4-13)  shows  that  summer  and  winter  cyclones 
have  similar  patterns  but  that  the  summer  patterns  are  naturally  less  intense. 


TABLE  4-3 

CHARACTERISTICS  OF  SUMMER  CYCLONES  OVER  NORTH  AMERICA, 
1955-1958  (DEPENDENT  SAMPLE) 


Forecast 
i nterva 1 , 
hr 

Observed 
northward 
d i sp lacement, 
deg  lat 

Observed 
eastward 
d i sp lacement, 
deg  lat 

Observed 
change  in 
central  pressure, 
mb 

Mean* 

Std  dev 

Mean  t 

Std  dev 

Meanf 

Std  dev 

12 

0.84 

2.04 

-5.50 

2.02 

-t  .08 

5.72 

24 

1.83 

5.58 

-6.4o 

5.69 

-1.69 

5.57 

*Nega+ive  values  represent  southward  displacement. 
tNegative  values  represent  eastward  displacement. 
{Negative  values  represent  deepening. 
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Fig.  4-15.  Mean  map  of  12-hr  pressure  change  for  summer 
sample).  Isallobars  are  labeled  in  millibars. 
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Fig.  4-17.  Mean  map  of  12-hr  500-mb  height  change  for  summer  cyclones  in  North  America,  1955-1958 
(dependent  sample).  Isallohypses  are  labeled  in  tens  of  feet. 


i 


n 


I 

I 

I 


k*IO 

Fig.  4-19.  Mean  map  of  12-hr  1000— 500-mb  thickness  change  for  summer  cyclones  in  North  America,  1955- 
1958  (dependent  sample).  Isallopachs  are  labeled  in  tens  of  feet. 


5.0  PREDICTION  EXPERIMENTS 


All  the  available  possible  predictors  of  each  predictand  must  be  reduced 
to  a  subset  to  be  used  in  the  screening-regression  prediction  experiments  because 
a  single  computer  run  cannot  use  all  3580  possible  predictors  (Table  3-2).  For 
each  predictand,  experimental  analysis  of  various  predictor  subsets  yielded  a  near¬ 
optimum  subset,  which  led  in  turn  to  a  "base -technique "equation  applicable  to 
the  entire  area  of  interest  and  to  all  cases  included  in  the  winter-cyclone  develop¬ 
mental  sample.*  Subsequent  analyses,  in  which  the  developmental  sample  was 
stratified  by  cyclone  position  at  forecast  time  and  additional  possible  predictors 
were  screened,  were  aimed  either  at  improving  the  results  of  the  base -technique 
equations  or  producing  equations  applicable  at  observation  times  other  than  those 
included  in  the  developmental  sample,  i.e.,  observation  times  at  which  only  sur¬ 
face  data  are  available. 

5.1  Base -technique  Prediction  Experiments 

In  the  initial  experiments,  predictors  of  only  types  A  and  B  (Table  3-2) 
were  considered.  The  multitude  of  possible  predictors  was  reduced  by  increasing 
the  spacing  between  gridpoints.  (Because  meteorological  observations  are  spatially 
redundant,  this  hypothetically  lost  very  little  of  the  predictive  information  contained 
in  the  complete  set  of  type  A  and  B  predictors.)  The  remaining  number  of  possible 
predictors  was  still  more  than  the  maximum  (170  to  175)  that  can  be  efficiently 
analyzed  by  screening  regression  (as  presently  programmed  for  the  IBM  7090), 
so  it  was  decided  to  conduct  two  base -technique  experiments.  Surface  (type  A) 
predictors  were  combined  with  700-mb  predictors  in  one  experiment  and  with 
500-mb  predictors  in  the  other.  The  combination  with  500-mb  predictors,  as 
listed  in  Table  5-1,  proved  slightly  better,  so  the  pertinent  equations  were 
defined  as  the  base  technique.  Results  of  the  base  technique  on  dependent  and 
independent  data  are  given  in  Section  6.0. 

♦Prediction  experiments  for  summer  cyclones  are  discussed  in  Section  5.6. 
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TABLE  5-1 

POSSIBLE  PREDICTORS  FOR  USE  IN 
BASE-TECHNIQUE  ANO 
STRATIFICATION  EXPERIMENTS 


Type 

Symbol* 

No.  usedt 

* 

1 

X 

1 

A 

P 

28 

6 P 

26 

2 

*10 

02 

37 

B 

H 

18 

_ 11 _ 

♦See  Table  5-2  for  dot inltions. 
tTotal  of  both  types  Is  168. 


TABLE  5-5 

POSSIBLE  PREDICTORS  FOR  USE  IN 
EXPERIMENTS  INVOLVING  SURFACE 
DATA  ONLY 


Type 

Symbol* 

No.  usedt 

♦ 

1 

). 

1 

A 

P 

» 

H- 

» 

♦See  Table  5-2  for  definitions* 
tTotal  is  17C. 


TABLE  5 m2 

POSSIBLE  PREDICTORS  FOR  USE  IN 
EXPERIMENTS  INVOLVING  DERIVED 
TERMS 


Type 

Symbol* 

No.  usedt 

6 

1 

A 

X 

1 

P 

12 

OP 

9 

2 

7 

CJ 

12 

h7 

8 

B 

*»7 

8 

H 

8 

4H 

7 

n 

11 

Aij 

10 

?T 

1C 

a*t 

9 

An 

1C 

ah7 

12 

C 

ah 

11 

A?T 

1C 

L 

9 

Uf|  V- 

r  < 

2 

Oujf  Oh- 

2 

e,  y,  IV i 

5 

u^f  V* 

1  4 

•Seo  Table  for  definitions. 
tTotal  of  ell  three  types  is  1  ;'l. 
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5.2  Geographic-stratification  Experiments 

Experiments  were  performed  to  determine  the  usefulness  of  equations 
derived  from  developmental  samples  stratified  according  to  position  of  the 
predictand  cyclone  at  forecast  time.  The  developmental  sample  was  stratified 
into  four  zones  (see  Fig.  2-1  and  Table  2-1).  The  possible  predictors  considered 
are  listed  in  Table  5-1.  Results  and  zone-by-zone  comparisons  with  the  base- 
technique  equations  are  presented  in  Section  6.0. 

5.3  Incorporation  of  Past-history  Predictors 

The  original  cyclone-data  tabulation  includes  a  6-hr  history  of  each  cyclones 
motion  (northward  and  eastward  components)  and  change  in  central  pressure.  To 
test  tiie  usefulness  of  these  three  possible  predictors,  they  were  added  to  the  base 
technique  and  the  prediction  experiment  was  repeated.  The  geographic- 
stratification  experiment  was  repeated  for  all  zones  and  all  predictands  because 
it  was  suspected  that  there  might  be  a  variation  from  zone  to  zone  in  the  accuracy 
with  which  the  past-history  predictors  could  be  estimated-for  example,  that  the 
position  and  central  pressure  of  cyclones  over  the  Atlantic  Ocean  and  Gulf  of 
Mexico  (in  the  southeastern  zone)  might  be  subject  to  greater  errors  than  those 
over  land  (in,  say,  the  northwestern  zone). 

The  possible  predictors  considered  in  these  experiments  are  listed  in 
Table  5-1,  with  the  3  past-history  predictors  included. 

5.4  Experiments  Using  Type  C  Predictors 

The  third  attempt  to  improve  the  base -technique  prediction  equations  was 
to  introduce  various  parameters  suggested  in  the  literature  as  being  important 
predictors  of  the  displacement  and  change  in  intensity  of  cyclones.  A  description 
of  this  class  of  predictors  is  contained  in  Section  3.4. 

The  required  reduction  in  the  number  of  possible  predictors  (Table  3-2) 
was  accomplished  by  computing  and  mapping  correlation  coefficients  of  all 
possible  predictors  with  each  of  the  predictands.  Predictors  then  selected  for 
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further  consideration  were  those  associated  with  extrema  of  the  mapped  corre¬ 
lation  fields  and  are  shown  in  Table  5-2. 

The  experiment  with  derived  terms  was  by  no  means  exhaustive,  but  it 
does  represent  an  attempt  to  incorporate  professional  knowledge  into  the 
statistical  analysis. 

5.5  Development  of  Prediction  Equations  from  Surface  Data  Only 

The  primary  objective  of  developing  prediction  equations  from  surface 
data  only  was  not  to  improve  the  results  obtained  by  the  base  technique  but  to 
increase  the  prediction  methods'  usefulness. 

Application  of  the  base  technique  is  limited.  Because  500-mb-height 
data  are  received  only  every  12  hr,  regression  equations  dependent  on  upper- 
air  data  can  be  applied  only  at  initial  times  of  0000  and  1200  GCT.  One 
approach  to  this  problem  is  the  use  of  500-mb  prognoses.*  Another  approach 
is  to  derive  regression  equations  that  use  only  surface  data  as  predictors. 

A  regression  analysis  that  excludes  upper-air  data,  if  successful,  would 
provide  a  useful  means  for  making  a  forecast  at  any  observation  time  for 
which  surface  data  are  available.  Accordingly,  two  experiments  were  per¬ 
formed.  All  cases  in  the  developmental  sample  were  analyzed.  The  possible 
predictors  for  these  experiments  are  listed  in  Table  5-3.  One  experiment 
included  the  past-history  predictors;  the  other  did  not.  Results  of  these  two 
experiments  are  presented  in  Section  6.0  for  dependent  and  independent  data. 

5.6  Prediction  Experiments  for  Summer  Cyclones: 

Four  experiments  were  conducted  to  develop  prediction  equations  for 
the  12-  and  24-hr  displacement  and  change  in  central  pressure  of  summer 
cyclones. 

The  first  experiment  employed  the  same  predictors  as  the  base-technique 
experiment  for  winter  cyclones  described  in  Section  5.1  (Table  5-1).  All 

♦The  authors  feel  that  this  approach,  not  yet  attempted,  should  be  attempted. 
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dependent  cases  (700)  in  the  entire  North  American  area  were  analyzed, 
and  the  resulting  equations  were  tested  on  independent  cases. 

The  second  experiment  included  the  past-history  predictors. 

Finally,  equations  were  derived  from  surface  data  only,  for  application 
when  upper-air  data  are  not  available.  As  in  the  previous  pair  of  experiments, 
the  first  prediction  experiment  did  not  include  the  past-history  predictors, 
and  the  second  one  did.  The  possible  predictors  are  those  listed  in  Table  5-3. 

Results  of  these  four  experiments  are  summarized  in  Section  6.0. 
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6.0  RESULTS 


The  first  series  of  results  discussed  (Section  6.1)  are  those  pertaining 
to  winter  cyclones.  Results  obtained  from  prediction  experiments  for  summer 
cyclones  are  presented  in  Section  6.2.  All  equations  developed  can  be  found 
in  Appendix  B. 

6.1  Winter  Cyclones 

The  sequence  of  experiments  for  winter  cyclones  was  first  to  establish 
base -technique  prediction  equations  for  the  north-south  (N)  and  east-west  (E) 
components  of  displacement  and  the  change  in  central  pressure  (D).  Equations 
were  derived  for  forecast  periods  of  12,  24,  and  36  hr.  The  entire  develop¬ 
mental  sample  was  used  to  derive  these  equations,  and  only  parameters 
routinely  available  at  0000  and  1200  GCT  were  considered  as  possible  pre¬ 
dictors.  Thus,  these  equations  are  applicable  to  all  North  American  cyclones 
regardless  of  location,  stage  of  development,  or  prevailing  synoptic  pattern. 
Development  of  these  equations  and  the  results  of  their  application  on  inde¬ 
pendent  data  are  given  in  Section  6.1.1. 

Sections  6.1.2  through  6.1.3  are  devoted  to  discussions  of  the  results  of 
experiments  aimed  at  improving  the  results  obtained  using  the  base-technique 
prediction  equations.  Improvements  were  sought  by 

(a)  development  of  equations  derived  from  developmental  samples 
stratified  according  to  the  location  of  the  cyclone  of  interest  at  forecast  time, 

(b)  introduction  of  the  displacement  and  change  in  central  pressure  of 
the  cyclone  during  the  6  hr  prior  to  forecast  time  as  possible  predictors,  and 

(c)  introduction  of  additional  derived  parameters  as  possible  predictors. 

Finally,  results  of  the  derivation  and  testing  of  equations  using  only  surface 
parameters  are  presented  in  Section  6.1.5. 
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6.1.1  Results  of  Uaingboth  Surface  and  Upper-air  Predictors 

Regression  equations  were  derived  for  the  prediction  of  the  12-,  24-, 
and  36-hr  displacement  and  change  in  central  pressure  for  North  American 
cyclones.  The  possible  predictors  are  listed  in  Tabic  5-1. 

The  results  of  the  screening-regression  analysis  are  summarized 
in  Tables  6-1  and  6-2.  Table  6-1  lists  the  predictors  in  the  order  of  their 
selection  by  the  screening  procedure  and  the  percentage  of  the  total  vari¬ 
ance  of  the  predictand  explained  by  each.  Table  6-2  summarizes  the  initial 
and  residual  standard  deviations  and  the  percent  reductions.  The  numbers 
accompanying  the  predictor  symbol  in  Table  6-1  refer  to  the  (k,f)-predictor 
locations  in  the  grid  system  shown  in  Fig.  3-1. 

From  Table  6-1,  it  can  be  seen  that  the  displacement  prcdictands  counted 
heavily  on  upper-air  information.  At  least  for  the  first  few  predictors  selected, 
the  sign  of  the  regression  coefficient  is  usually  consistent  with  qualitative 
reasoning.  For  example,  in  each  case,  the  second  predictor  selected  for 
the  prediction  of  the  cyclones'  northward  component  of  displacement  was  a 
500-mb  height  from  3  to  5  grid  intervals  cast  of  the  storm.  The  regression 
coefficient  associated  with  this  predictor  is  positive  in  all  three  northward- 
displacement  equations:  thus,  a  large-amplitude  ridge  east  of  the  cyclone 
would  contribute  to  the  prediction  of  a  sizable  northward  displacement. 

It  is  interesting  to  note  that  the  longitude,  X,  of  the  cyclone  center  at 
forecast  time  is  the  first  predictor  selected  for  each  of  the  three  northward- 
displacement  equations.  The  sign  of  the  regression  coefficient  (negative) 
indicates  that  this  is  likely  a  reflection  of  the  variation  of  the  climatology 
of  storm  tracks  within  the  prediction  region,  i.e.,  a  tendency  for  a  south¬ 
ward  component  of  displacement  in  the  west  and  for  a  northward  component 
in  the  east. 

The  first  two  predictors  selected  for  the  three  eastward-displacement 
equations  follow  the  same  pattern.  The  first  predictor  selected  is  a  500-mb 
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TABLE  6-1 

PREDICTORS  SELECTED  RY  SCREEN  INC  REGRESSION* 
FOR  WINTER  CYCLONES  OVER  NORTH  AMERICA 
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height  north-northeast  of  the  surface  cyclone  center  and  the  second,  a  500-mb  height 
south-southwest  of  the  center.  The  regression  coefficients  associated  with  the  two 
predictors  are  approximately  equal  in  magnitude  and  of  opposite  sign— in  effect, 
a  measure  of  the  strength  of  the  eastward  "steering"  component. 

The  central  pressure's  lead  predictors,  AP(9,5)  for  12  hr  and  AZ(9,.r>)  for  24  hr, 
correlate  12-hr  pressure  and  height  falls  to  the  immediate  west  of  the  cyclone  with 
deepening.  The  first  two  predictors  selected  for  the  36-hr  change-in-central-pressure 
equations  and  their  coefficients  (opposite  in  sign)  indicate  a  correlation  between  the 
strength  of  the  westerlies  and  deepening. 

Table  6-3  contains  the  results  of  application  of  the  base-technique  equations 
to  213  cases  not  used  in  their  development.  A  comparison  of  independent-data  rms 
errors  with  the  residual  standard  deviationsof  the  developmental  sample  indicates  a 
remarkable  degree  of  stability  in  the  equations. 

6.1.2  Results  of  the  Geographic  Stratification  of  the  Developmental  Sample 

Geographic  variations  in  cyclone  climatology  and  in  the  accuracy  with  which 
predictors  may  be  specified  suggest  that  subdividing  the  developmental  sample 
according  to  the  location  of  the  predictand  cyclone  and  rederiving  equations  appli¬ 
cable  to  zones  within  the  forecast  area  may  lead  to  improved  results.  Nevertheless, 
unless  very  large  samples  are  available,  errors  in  estimating  regression  eoefii- 
cients  employing  subsets  of  the  entire  developmental  samples  may  counteract  the 
improvement. 

To  investigate  this  possibility  for  improving  the  base -technique  equations, 
separate  regression  equations  were  derived  for  cyclones  by  zone  (Fig.  2-1).  The 
predictands  and  possible  predictors  were  exactly  the  same  as  in  the  previous 
experiment,  which  considered  all  cyclones  in  the  entire  forecast  area. 

The  results  of  the  screening-regression  analysis  are  summarized  in  Tables 
6-4  and  6-5.  Table  6-4  lists  the  predictors  in  the  order  of  their  selection  bv  the 
screening  and  the  percentage  of  the  total  variance  of  the  predictand  explained  by 
each.  Table  6-5  summarizes,  by  zone,  the  initial  and  residual  standard  deviations 
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TAfJI  E  6-4 
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FOR  WINTER  CYCI ONES  OVER  NORTH  AMERICA 
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H(  7, 

9) 

3.6 
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- 

- 

*>(  7, 

5) 

3.1 

Total 

53.6 

59.1 
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9) 
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3.5 

az(  9, 

9) 

4.6 

2U 
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7(15,  7) 

n 

7(  7.  5) 
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P(  5, 

3) 

3.1 

6th 

H(  7,  9) 

E9 

7(15,11 ) 

1.6 
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1) 

1.8 

7th 

- 

■ 

az(  5,  3) 

1.7 

aZ(ll. 

7) 

0.9 

8th 

- 

- 

• 

aP(  9, 

3) 

0.8 

9th 

- 

■ 

- 

- 

P(10, 

5) 

1.3 

10th 

- 

■i 

- 

- 

*11. 

7) 

2.4 

i  Total 

56.6 

mmmmsm 

55.4  | 
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29.0 

9, 

3) 
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Z(13,  7) 

28.1 
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D 

11.0 

1/ 

3rd 
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4.3 
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9) 

1 

x> 

4th 
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Z(  9, 
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H 

>th 
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3.0 
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5) 

D 
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- 

- 
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1.8 

- 

Wm 

|  Total 

1  53.0 

62.6 

43.9 
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•Including  upper-sir  predictors 
tsijoifieant  at  the  5'<'  level. 


•Inelulinj  upper-#ir  predictors, 
t^ijnifieant  *t  the  5';  level. 


and  the  percent  reductions.  The  numbers  accompanying  the  predictor  symbols  in  Table 
6-4  refer  to  the  (k,l)-predictor  locations  in  the  grid  system  shown  in  Fig.  3-1. 

Table  6-6  contains  the  results  of  the  application  of  the  equations,  by  zone, 
to  the  independent  sample  of  213  cases.  Also  listed  in  Table  6-6  are  the  rms 
errors,  by  zone,  resulting  from  the  application  of  the  base -technique  equations. 

Although  11  of  the  36  stratified  equations  resulted  in  lower  rms  errors  than 
those  associated  with  the  same  forecasts  using  the  base -technique  equations,  only 
the  northwestern-zone  equation  for  the  36-hr  northward  displacement  was  signif¬ 
icantly  better  (at  less  than  the  5%  level)  than  the  base -technique  equation.  The 
significance  test  employed  is  Student's  t-test  for  paired  comparisons,  described 
by  Wadsworth  and  Bryan  ( 15] . 

Accordingly,  the  base  technique  was  redefined  for  predictions  of  the  36-hr 
northward  displacement  in  the  northwestern  zone. 

6.1.3  Results  of  Incorporating  Past  History  as  Possible  Predictors 

The  next  experiments  were  carried  out  to  see  if  improvements  could  be 
realized  by  using  the  past  history  as  predictors.  Specifically,  three  new  possible 
predictors  were  considered.  They  are  the  north-south  and  east-west  components 
of  a  cyclone's  displacement  during  the  6-hr  period  just  prior  to  forecast  time 
and  the  change  in  central  pressure  of  the  cyclone  during  the  same  time.  The  sole 
difference  between  these  experiments  and  the  previous  series  was  the  consideration 
of  the  three  additional  possible  predictors.  Both  the  stratified  and  unstratified 
samples  were  analyzed  because  of  the  possibility  that  errors  in  estimating  the 
past  track  and  change  in  central  pressure  would  vary  from  zone  to  zone. 

The  results  of  the  screening-regression  analysis  are  summarized  in  Tables 
6-7  and  6-8.  Table  6-7  lists  the  selected  predictors,  in  order  of  selection  by 
the  screening  procedure,  and  the  percentage  of  the  total  variance  explained  by 
each.  Table  6-8  summarizes  the  initial  and  residual  standard  deviations  and  the 
percent  reductions. 
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7ABI  E  6-  ; 

PREDICTORS  SELECTED  BY  SCREEN  I  tIG  REGRESSION 
FOR  WINTER  CYC  I  ONES  OVER  NORTH  AMERICA* 
(PAST  HISTORY  INCtUDED) 

(a)  Nor thwuulorn  ?ono 


Forecast 
interval , 
hr 

Order  ot 
select  ion 

N 

E 

6 

Predictor 

1 

red 

r  red ictor 

I  ' 

red 

l  red ictor 

HI 

1st 

*\0 

m 

22.1 

07 (  5,  7) 

5.2 

2nd 

0vH(1 3,  0) 

7(11,  9) 

10.0 

P(1C,  5) 

10.3 

12 

3rd 

7(13,  5) 

i 

7(13,  1) 

8.2 

V(  9,  5) 

9.o 

4th 

7(  9,  7) 

4.1 

- 

- 

*(11,  5) 

11.3 

Oth 

> 

Vv 

- 

- 

- 

- 

j  Total 

v8.i 

■It.  3 

^  1 

1  St 

«.  0 

26.3 

7(11,  9) 

24.4 

'(  5,  7) 

12.2 

2nd 

'-VI(13,  9) 

7.9 

’(11,  1) 

15.9 

P(1o,  -j) 

15.5 

3rd 

X 

,J*7 

V(  9,n) 

3.7 

V(  9,  5) 

"•‘t 

24 

4th 

7(19,  3) 

5*5 

1(13,  7) 

2.7 

F(11,  5) 

6.5 

Oth 

V(13,  9) 

3.4 

’(  9,  3) 

5.7 

■^e 

•«  •  r 

6th 

7(  9,  7) 

3.2 

- 

- 

- 

- 

7th 

7(13,  9) 

4.9 

- 

- 

- 

- 

|  Tcral 

59-9 

47.6  | 

1st 

H(!3,  5) 

1 7.3 

7(11,  9) 

25.0 

p(lt,  5) 

12.7 

56 

2nd 

"(  7,  7) 

7(11,  1) 

16.5 

■•’(  9,  7) 

1 3.2 

3rd 

- 

- 

- 

- 

7(15,  0 

6.u 

|  Total 

33.0 

’*1.5 

51.9  | 

•Including  upper-air  predictors. 
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Table  6-9  shows  the  rms  errors,  by  zone,  found  in  applying  the  two  sets 
(stratified  and  unstratified)  of  derived  equations.  For  comparison,  the  rms  errors 
resulting  from  the  application  of  the  base -technique  equations  to  the  same  sample 
are  included. 

Student's  t-test  for  paired  comparisons  indicated  that  equations  that  include 
past  history  as  predictors  yield  improvements  in  two  instances  that  are  statistically 
significant  at  less  than  the  5%  level.  One  is  the  application  to  the  southwestern  zone 
of  the  12-hr  eastward-displacement  equation  derived  from  the  unstratified  develop¬ 
mental  sample.  The  other  is  the  application  to  the  northeastern  zone  of  the  24-hr 
northward-displacement  equation  derived  from  the  stratified  sample. 

6.1.4  Results  of  Incorporating  Additional  Derived  Predictors 

After  the  necessary  reduction  of  possible  derived  predictors  was  made  by 
the  method  described  in  Section  5.0,  two  screening-regression  experiments  were 
performed.  Results  of  the  test  on  independent  data  were  compared  with  those  of 
the  modified  base  technique. 

The  two  experiments  employing  additional  derived  predictors  were  limited 
to  consideration  of  the  three  24-hr  predictands,  and  equations  were  derived  using 
only  the  unstratified  developmental  sample.  The  first  screening-regression 
experiment  considered  the  possible  predictors  listed  in  Table  5-2.  The  second 
experiment  considered,  in  addition,  the  three  past-history  predictors. 

The  results  of  the  screening-regression  analysis  are  summarized  in  Tables 
6-10  and  6-11.  Table  6-10  lists  the  predictors  in  the  order  of  their  selection  by 
the  screening  procedure  and  the  percentage  of  the  total  variance  of  the  predlctand 
explained  by  each.  Table  6-11  contains  a  summary  of  the  initial  and  residual 
standard  deviations  and  the  percent  reductions. 

Table  6-12  summarizes  the  results,  by  zone,  of  the  application  of  the  equations 
to  the  independent  data  sample  of  213  cases.  Also  included  are  the  results  obtainsd 
by  the  modified  base  technique. 
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TABLE  6-10 

PREDICTORS  SELECTED  BY  SCREENING  REGRESSION 
FOR  WINTER  CYCLONES  OVER  NORTH  AMERICA 
(INCLUOING  DERIVED  PREDICTORS) 


Forecast 
interval , 
hr 

Order  of 
selection 

N 

E 

D 

Predictor 

% 

red 

Predictor 

% 

red 

Predictor 

(  ' 

<( 

red 

1st 

v7 

33.6 

u7 

29.3 

1  v  1 

19.8 

2nd 

X 

8.2 

2(12,10) 

10.1 

X 

7.9 

3rd 

P<16,  6) 

4.1 

A„  (9,7) 

4.3 

a2(  9,  8) 

5.6 

lith 

*H7 

3.3 

a2(  6,  5) 

2.9 

P(10,  5) 

4.8 

5th 

2.5 

2(11,  D 

1.9 

az(to,  3) 

3.4 

24* 

6th 

p(io.  6) 

2.8 

v7 

2.0 

ah  (6,4) 

3.1 

7th 

<t(8.  8) 

1.1 

atT('0,9) 

1.2 

a?t(9,io) 

1.4 

8th 

Z(l4,  4) 

1.1 

p(  5,  5) 

1.0 

*h  (9,  7) 

1.4 

9th 

«  (6.  8) 

1.9 

p(i  7,  3) 

1.0 

aP(9,  6) 

2.9 

10th 

H?(12,2) 

1.2 

act(13,4) 

v7 

1.2 

11th 

- 

- 

- 

m 

5P(1C,  3) 

1.2 

|  Total 

59.8 

54.6 

»■’  1 

1st 

a* 

39.7 

OA 

30.5 

1  v  1 

19.8 

2nd 

7(14,  4) 

7.4 

u 

8.9 

X 

7.9 

3rd 

n  i,  7) 

6.3 

2(12,10) 

6.1 

V(  9,  8) 

4th 

X 

2.8 

v7 

3.4 

P(10,  5) 

K| 

5th 

CT(S»  8) 

1.7 

h  (6,  8) 

2.2 

02(10,  3) 

111 

6th 

P(10,  6) 

1.3 

5CT(,0*9) 

1.1 

3.1 

7th 

H7(12,2) 

1.3 

P(13,  1) 

1.0 

2.1 

24t 

8th 

*7 

1.2 

r— 

ON 

X 

< 

0.8 

4?T(9,10) 

1.2 

9th 

l»l 

0.9 

H7(12.2) 

0.8 

*„  (9,  7) 

1.4 

10th 

CT(14,12) 

0.7 

v7 

■a 

3P(  9,  6) 

1.6 

11th 

- 

- 

P(  5,  9) 

1 

* 

v7 

1.1 

12th 

- 

- 

# 

y7 

0.8 

DP(10,3) 

1.0 

13th 

- 

- 

^7(6,9) 

0.9 

- 

- 

Total 

63.3 

58.3 

53.0  | 

•excluding  past  history, 
tlncluding  past  history. 
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An  examination  of  Table  6-12  reveals  improvement  over  the  base  technique  in 
five  of  the  12  comparisons.  However,  only  application  to  the  northwestern  zone  of 
the  derived-term  (without  past  history)  equation  for  24-hour  northward  displacement 
yielded  an  improvement  statistically  significant  at  less  than  the  5%  level.  The  base 
technique  was  modified  to  include  this  equation. 

6.1.5  Results  of  the  Application  of  Prediction  Equations  Using  Surface  Data  Only 

To  make  the  prediction  technique  applicable  to  those  synoptic  observation  times 
for  which  upper-air  observations  are  not  available,  equations  were  derived  from  sets 
of  possible  predictors  containing  only  surface  data.  Two  screening-regression 
experiments  excluding  possible  upper-air  predictors  were  carried  out.  The  exper¬ 
iments  were  identical  except  that  three  additional  possible  predictors  (past  history) 
were  included  in  the  second  experiment. 

In  these  experiments,  prediction  equations  were  derived  for  a  cyclone's  12- 
and  24-hr  displacements  and  changes  in  central  pressure.  The  possible  predictors 
common  to  both  experiments  are  given  in  Table  5-3. 

Results  of  the  screening-regression  analysis  are  summarized  in  Tables 
6-13  and  6-14.  Table  6-13  lists  the  predictors  in  the  order  of  their  selection  by 
the  screening  procedure  and  the  percentage  of  the  total  variance  of  the  predictand 
explained  by  each.  Table  6-14  summarizes  the  initial  and  residual  standard 
deviations  and  the  percent  reductions. 

Table  6-15  summarizes  the  results  of  testing  the  equations  derived  from 
surface  predictors  on  the  independent  data  sample.  The  results  obtained  from 
the  base -technique  equations  are  also  included. 

Although  the  overall  performance  of  the  base  technique  is  slightly  superior 
to  that  of  the  equations  derived  only  from  surface  predictors,  it  appears  that  useful 
results  can  be  obtained  from  the  surface -predictor  equations  when  upper-air  data 
are  not  available.  Note  from  Table  6-15  that  the  12-hr  change-in-central-pressure 
equation  based  only  on  surface  data  and  past  history  yielded  a  considerably  smaller 
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TABLE  6-1} 

PREOICTORS  SELECTED  BY  SCREENING  REGRESSION* 
FOR  WINTER  CYCLONES  OVER  NORTH  AMERICA 

(a)  Excluding  past  history 


Forecast 
interval , 
hr 

Order  of 
selection 

N 

- i 

6 

Predictor 

% 

red 

Predictor 

% 

red 

Predictor 

% 

red 

1st 

X 

22.6 

P(  1 3 , 1 1 ) 

1 4 .6 

4P(  9,  5) 

MH 

2nd 

P(15, 

7) 

8.7 

i>P(  7.  5) 

7.4 

X 

3rd 

P(10, 

6) 

7.0 

P(11,  0 

4.6 

p(io,  5) 

7 

4th 

aP(  9, 

7) 

2.8 

P(  6,  6) 

2.5 

P(n,  5) 

5th 

op(  7, 

3) 

2.7 

P03,  7) 

2.5 

4P(  5,  3) 

6th 

P(ll, 

7) 

2.3 

P(17,  3) 

3.0 

6P(  10,  4) 

12 

7th 

P(  6, 

2) 

2.3 

P(  9,  3) 

1-9 

0P(10,  8) 

2.1 

8th 

P(13» 

3) 

1.6 

4P(10,  4) 

2.7 

P(17»  5) 

1.2 

9th 

0P(  10, 

2) 

1.8 

p(  7,  O 

1.8 

P(  6,  6) 

0.8 

10th 

0P(11, 

9) 

1.0 

P(13,  5) 

1-7 

*>(  8,  6) 

1.1 

11th 

p(  5, 

9) 

1.1 

P{14,  4) 

0.9 

- 

- 

12th 

P(15, 

D 

0.7 

*P(  8,  6) 

0.8 

- 

- 

nth 

- 

- 

4P(  3,  5) 

1.0 

- 

- 

Total 

54.6 

45.4 

42.9 

1st 

X 

22.3 

P(15,11) 

16.9 

X 

15.4 

2nd 

P(17, 

7) 

7.8 

P(13,  1) 

7-5 

0P(  9,  3) 

7.0 

3rd 

P(10, 

6) 

6.0 

P(13,  7) 

3.1 

0P(  5,  3) 

4.5 

4th 

P(l4, 

4) 

3.6 

X 

3.1 

P(10,  5) 

5.5 

5th 

P(  7, 

1) 

2.5 

♦ 

3.7 

0P(10,  6) 

4.5 

24 

6th 

0P(  7, 

1) 

2.4 

AP(  7,  5) 

1.8 

P(12,  6) 

4.8 

7th 

6P(  9, 

7) 

1.6 

P(  5,  7) 

2.1 

p07»  3) 

2.5 

8th 

p(ll, 

7) 

2.2 

P(10,  2) 

2.5 

ap(io,  4) 

1.1 

9th 

0P(11, 

9) 

1.0 

flP(11»  5) 

1.3 

p(ii,  1) 

1.8 

10th 

P(  1» 

7) 

1.1 

P(17,  3) 

1.4 

- 

- 

11th 

P(l5, 

7) 

1.1 

P(  9/  7) 

0.9 

- 

- 

Total 

51.6 

45.8 

45.1  ] 

♦Excluding  upper-air  pradictors. 
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(b)  Including  past  history 


Foracast 

intarval, 

hr 

Ordar  of 
salactlon 

N 

E 

6 

Pradictor 

% 

rad 

Pradictor 

% 

rad 

Pradictor 

% 

rad 

i.t 

A# 

51.4 

AX 

e§hs 

9. 5) 

■fli 

2nd 

X 

).o 

p(i),ii) 

X 

El 

irt 

p07,  7) 

2.0 

P(10,  2) 

1 

p(10,  5) 

4.7 

4th 

p(ii,  7) 

2.1 

P(12,  6) 

111 

P(11,  5) 

6.2 

5th 

P04,  6) 

1.5 

X 

1.7 

9P(  5,  5) 

5.5 

12 

6th 

♦ 

1.2 

P(15,  5) 

1.5 

xOc.,  ** ) 

5.0 

7th 

.P05,  5) 

0.6 

9P(  8,  6) 

1.0 

At 

5.7 

8th 

5P(  9,  )> 

0.7 

* 

1.1 

9»'c 

9th 

- 

- 

p(  5,  7) 

1.5 

- 

■ 

10th 

- 

- 

5P{1I,  5) 

1.0 

- 

11th 

- 

- 

p(10,  6) 

1.5 

- 

H 

Total 

62.) 

52.5 

1*5.4  | 

1st 

A# 

59.7 

AX 

50.5 

X 

15.4 

2nd 

X 

■*.5 

pO),h) 

6.5 

X(  9,  5) 

7.0 

yd 

p(h,  ;) 

4.0 

pO  5,  1) 

2.4 

'H  5,  5) 

4.5 

4th 

p(n,  7) 

5.9 

X 

2.1 

pO0,  5) 

5.5 

5th 

V(15.  5) 

1.2 

0 

2.2 

4.  ,• 

24 

6th 

♦ 

1.2 

P(15.7) 

2.8 

P(11,  5) 

5.5 

7th 

p04,  •>) 

0.8 

p(  5.  )) 

1.2 

V(10,  5) 

4.5 

8th 

)P(n,  9) 

0.6 

p(W,  5) 

1.0 

A# 

1.7 

9th 

ap(H>,  6) 

u.7 

'(  9,  5) 

1.1 

p(n.  7) 

1.1 

10th 

5P(tl,  )) 

O.t 

V'O),  9) 

0.9 

p(l 7.  5) 

1.5 

11th 

- 

- 

T(  8,  6) 

0.8 

- 

- 

12th 

- 

- 

P(  6,  8) 

1.2 

- 

* 

Total 

57.) 

52.5 

*»7.o  | 

TABLE  6-l4 

RESULTS  or  SCREENING  REGRESSION*  ON  WINTER  CYCLONES  OVER  NORTH  AFRICA, 
195V 1959  ( DEPENDENT  SAfcTlE) 


7ona 

Foracast 

intarval, 

No.  Of 
pradictors 

Std  dav 

Rasidual  std 

dav 

^  raduction 

KB 

E 

0 

N 

E 

6 

N 

E 

6 

N 

£ 

6 

1? 

12 

15 

10 

2.69 

2.42 

5.55 

1.81 

1.56 

4.19 

54.6 

45.4 

42.9 

A 1 1 1 

24 

11 

11 

9 

4.69 

4.46 

9.54 

5.27 

5.28 

7.07 

51.6 

45.8 

45.1 

12 

8 

11 

6 

2.69 

2.42 

wvm 

1.65 

1.67 

4.17 

62.5 

52.5 

45.4 

A 1  It 

mm 

24 

10 

12 

10 

4.69 

4.46 

WSM 

5.07 

5.07 

6.94 

57.5 

52.5 

47.0 

•Excluding  uppar-air  pradictor*. 
tExcludlng  past  history. 
{Including  past  history. 
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rms  error  in  the  northeastern  zone  than  the  base  technique  did.  Student's  t-test 
for  paired  comparisons  indicates  that  the  improvement  is  significant  at  less 'than 
the  5%  level.  The  inclusion  of  this  equation  represents  the  final  modification  made 
to  the  base  technique  in  this  study. 

.  .  6.1.6  Operational  Tests 

To  assess  the  accuracy  of  the  regression  technique,  the  unimproved  base- 
technique  equations  were  applied  to  all  cyclones  identified  in  the  area  defined  by 
Fig.  2-1  between  November  1962  and  March  1963,  and  the  results  were  compared 
with  NWAC  24-  and  36-hr  prognoses  received  by  facsimile.  The  vector-position 
rms  errors  and  the  central-pressure  rms  errors  are  shown  in  Table  6-16.  Included 
also  are  the  results  of  the  stratified-sample  equations.  NWAC  prognoses  were  available 
for  57  of  the  eighty-one  24-hr  forecasts*  and  30  of  the  sixty-five  36-hr  forecasts. 

The  table  indicates  that,  for  this  sample,  the  base  technique  continued  to 
perform  well  in  comparison  with  stratification.  The  36-hr  reduced-sample 
forecast  of  central  pressure  is  the  only  important  exception.  The  comparisons 
with  the  NWAC  prognoses  Indicate  that  the  base  technique  was  better  for  24-hr 
position  and  36-hr  pressure,  but  the  NWAC  prognoses  were  better  for  24-hr 
pressure  and  36-hr  position. 

6.2  Simmer  Cyclones 

Prediction  experiments  similar  to,  but  fewer  than,  the  winter  series  were 
conducted  to  develop  prediction  equations  for  the  behavior  of  summer  cyclones 
over  North  America.  The  number  of  cases  available  for  development  of  the 
summer  equations  and  the  number  of  independent  cases  for  testing  the  results  are 
given  in  Table  2-1. 

Two  sets  of  prediction  equations  were  derived  for  application  at  observation 
times  for  which  both  surface  and  upper -air  observations  are  available;  their 

♦In  some  cases  the  cyclone  under  consideration  was  off  the  edge  of  the  NWAC 
facsimile  map. 
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TABLE  6*17 

PREDICTORS  SELECTEO  BT  SCREENING  REGRESSION* 
FOR  SUMCR  CYCLONES  OVER  NORTH  AICRICA 

(a)  Excluding  poit  history 


Forecast 

interval, 

hr 

Ordar  of 
»• tact  Ion 

N 

E 

D 

Pradlctor 

% 

rad 

Pradlctor 

m 

Pradietor 

% 

rad 

lit 

0P(  9,  7) 

9.3 

7(11,  9) 

16.4 

*>(  9,  5) 

2nd 

p<15,  7) 

7.0 

7(  9,  3) 

19.4 

P(10,  5) 

Jrd 

P(it.  7) 

8.0 

7(11,  7) 

5.3 

P(11,  5) 

5.1 

4th 

aP(  7,  5) 

3.2 

H(1S,  5) 

2.1 

«P(  7,  5) 

5.5 

5th 

flPOi,  9) 

2.1 

7(  3,  7) 

1.7 

0P(10,  5) 

2.7 

6th 

X 

2.1 

P(11,  1) 

1.6 

5P(  7,  5) 

1.4 

7th 

p(  9,  1) 

2.1 

aP(ii,  7) 

o.e 

- 

12 

8th 

7(  9,  7) 

1.3 

aP(  9,  5) 

1.4 

9th 

7(13,  5) 

13.3 

02(11,  3) 

0.7 

10th 

H(13,  3) 

1.5 

V(  7,  5) 

1.0 

11th 

P(15*  9) 

1.3 

- 

- 

12th 

M<  9,  5) 

1.1 

- 

- 

15th 

P(  5,  1) 

- 

- 

l4th 

*><10,  5) 

- 

- 

15th 

0P(11,  3) 

1.1 

- 

• 

- 

- 

Total 

55.2 

50.4 

26.0  | 

lit 

P(15,  7) 

jf§nn 

2(11,  9) 

18.5 

p(io,  5) 

7.3 

2nd 

p(",  7) 

1 

7(  9,  3) 

22.9 

0P(10,  5) 

10.4 

5r  d 

0P(  9,  7) 

7(11,  7) 

3.1 

P(11,  5) 

7.4 

4th 

«P<  7,  5) 

7(  3,  7) 

1.7 

*>(  7,  3) 

2.8 

>th 

X 

EH 

h(13,  5) 

2.5 

6th 

p<  9,  1) 

2.8 

P(15,  3) 

1.6 

7th 

N(  9,  7) 

2.0 

7(11,  5) 

0.9 

24 

8th 

7(13,  5) 

15.0 

07(11,  7) 

1.2 

9th 

7(15,11) 

1.8 

v(  7,  5) 

1.8 

10th 

7(  7,  1) 

1.7 

- 

11th 

«(  5,  1) 

1.1 

- 

12th 

aP(lo,  5) 

0.8 

- 

15th 

n(13,  3) 

0.7 

- 

14th 

5P(11,  5) 

1.0 

- 

|  Total 

58.0 

54.2 

27.9  | 

•Including  uppar.air  pradlctors. 
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I 

I 

1 

I 

I 

I 

I 

1 

I 
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(b)  Including  past  history 


Forecast 
intorva 1 , 
hr 

Or dor  of 
seloct ion 

N 

E 

0 

Predictor 

% 

red 

Predictor 

% 

red 

Pred ictor 

% 

red 

1st 

A  <t> 

29.2 

AX 

29.8 

AP(  9,  9) 

7.9 

2nd 

P(15,  7) 

3.1* 

7(11,  9) 

9.2 

AX 

4.1 

5rd 

P(11,  7) 

4.1 

7(  9,  5) 

8.9 

P(10,  9) 

3.9 

4th 

AP(  9,  7) 

2.9 

7(11,  7) 

2.6 

f’(11,  9) 

5.1 

91  h 

AX 

1 .4 

Af(  9,  9) 

1.9 

*p(  7,  3) 

4.5 

6th 

AP(15,  9) 

1.1 

Af'03,  9) 

1.6 

AF(10,  9) 

1.8 

— 

/th 

X 

0.8 

AF’(  9,  3) 

0.7 

- 

- 

if 

8th 

P(  9,  D 

1 .2 

*7(  3,  7) 

0.7 

- 

- 

<Hh 

P(13,  » 

0.9 

H(13,  9) 

1.2 

- 

- 

10th 

At'(11,  3) 

0.9 

- 

- 

- 

- 

11th 

h(  9,  7) 

o.9 

- 

- 

- 

- 

M 

12th 

?<15,  7) 

7-9 

- 

- 

- 

- 

|  Total 

94.3 

92.2 

26.7  | 

1st 

A* 

22.6 

AX 

23.4 

P(10,  9) 

7.3 

2nd 

p(19,  7) 

9.2 

7(11,  9) 

7.8 

AP(l  0,  9) 

io.4 

5rd 

P(",  7) 

6  •!> 

7(  9,  5) 

15.7 

P(11,  9) 

7*4 

lith 

7(  9,  7) 

2.1  • 

7(11,  7) 

1.7 

*T(  7,  5) 

2.8 

9th 

?05,  9) 

15.6 

AT(  9,  5) 

1.8 

P(11,  9) 

1.1 

24 

6th 

/0%11) 

1.8 

AP(13,  9) 

2.0 

*p(1 9,  7) 

0.8 

,’th 

♦ 

2.0 

*7(11,  7) 

1.1 

A0 

0.7 

8th 

P(  9,  1) 

1 .0 

V(  7,  9) 

1.4 

H(  7,  7) 

0.6 

9th 

H(  9,  9) 

1.2 

- 

- 

H(  11,  1) 

1.4 

10th 

AP(1I,  7) 

1.1 

- 

- 

v(  9,  9) 

1.5 

_ 

11th 

- 

- 

- 

- 

V(  9,  D 

3.5 

j  Total 

!  97*1 

92.9 

57.3 
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derivation  and  testing  are  discussed  in  Section  6.2.1.  Two  other  sets  of  prediction 
equations  were  derived  for  application  at  observation  times  for  which  only  surface 
data  are  available;  their  results  are  discussed  in  Section  6.2.2. 

6.2.1  Results  of  Using  Both  Surface  and  Upper-air  Predictors 

The  first  equations  derived  were  to  be  applicable  to  all  cyclones  in  the  entire 
North  American  prediction  area  (Fig.  2-1)  for  0000  and  1200  GCT  when  both  upper- 
air  and  surface  data  are  available.  To  make  the  equations  applicable  to  newly  formed 
cyclones,  the  motion  and  change  in  central  pressure  of  the  cyclone  just  prior  to 
forecast  time  were  not  considered  as  possible  predictors.  The  predictands  for  all 
summer-cyclone  prediction  experiments  were  the  12-  and  24-hr  displacement 
components  and  the  12-  and  24-hr  change  in  central  pressure.  The  possible 
predictors  for  this  first  experiment  are  listed  in  Table  5-1.  The  second  summer- 
cyclone  experiment  using  both  surface  and  upper-air  predictors  wsb  the  same  as 
the  first  except  that  the  three  past-history  predictors  (Table  3-2)  were  included. 

The  results  of  these  two  screening-regression  analyses  are  summarized  in 
Tables  6-17  and  6-18.  Table  6-17  lists  the  predictors  in  the  order  of  their 
selection  and  the  percentage  of  the  variance  of  the  predictand  explained  by  each. 

Table  6-18  summarizes  the  results  of  application  of  the  derived  equations  to  the 
developmental  sample. 

Table  6-19  lists  the  rms  errors  resulting  from  the  application  of  the 
equations  to  the  independent  data  sample  (169  cases).  The  rms  errors  associated 
with  the  equations  using  past-history  predictors  are  only  slightly  lower  than  those 
obtained  from  application  of  the  equations  not  using  past-history  predictors. 

6.2.2  Results  of  Using  Surface  Predictors  Only 

To  make  the  technique  applicable  to  observation  times  for  which  only 
surface  data  are  available,  two  sets  of  equations  were  derived  with  only  surface 
parameters  considered  for  screening.  The  set  of  possible  predictors  for  the  first 
experiment  is  listed  in  Table  5-3.  The  second  experiment  using  only  surface 
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TABLE  6-20 

PREDICTORS  SELECTED  BY  SCREENING  REGRESSION* 
FOR  SUMER  CYCLONES  OVER  NORTH  AMERICA 


(a)  Excluding  past  history 


Forecast 
interval , 
hr 

Order  of 
selection 

N 

E 

6 

Predictor 

% 

red 

Predictor 

% 

red 

Predictor 

% 

red 

1st 

AP(10,  6) 

13.0 

P(12,  8) 

Bii* 

0P(  9,  5) 

7-9 

2nd 

P(15,  7) 

7-2 

P(11,  D 

P(  1,  5) 

5.3 

3rd 

P(n,  7) 

4.2 

AP(  8,  6) 

BBlK 

P(10,  8) 

4th 

AP02,  8) 

3.5 

tf’dl,  5) 

P(10,  5) 

5th 

AP(  9,  7) 

2.7 

P(10,  5) 

P(ll,  5) 

12 

6th 

X 

2.0 

P(12,  6) 

4.o 

*>(  8,  4) 

^^9 

7th 

p(  9,  D 

1.9 

- 

- 

0P(1 5,  3) 

1.8 

8th 

&(  8,  4) 

1.6 

- 

- 

6P(10,  4) 

1.4 

9th 

P(11,  3) 

1.5 

- 

- 

- 

- 

10th 

p(15*11 ) 

1.6 

- 

- 

- 

- 

Total 

39.0 

33.9 

29.5 

1st 

P(15,  7) 

10.3 

P(13,  9) 

13.4 

P(10,  5) 

7.3 

2nd 

P(10,  6) 

10.6 

P(12»  2) 

7-4 

3P(10,  5) 

10.4 

3rd 

X 

3.9 

*>(  8,  6) 

3.7 

P(11»  5) 

7.4 

4th 

0P(  9,  7) 

2.4 

OP(1l,  7) 

3.1 

7,  5) 

2.8 

5th 

ar(  7#  5) 

3.4 

P(13,  5) 

3.3 

- 

- 

6th 

P(»5,11) 

1.9 

X 

1.3 

- 

- 

24 

7th 

p(l4,  4) 

1.8 

p(  1,  7) 

1.3 

. 

- 

8th 

P(12,  6) 

2.0 

P(  5,  3) 

1.1 

- 

- 

9th 

P(  9,  D 

2.1 

4P(11,  5) 

1.3 

- 

- 

10th 

P(  9,  9) 

1.8 

P(10,  4) 

2.1 

- 

- 

11th 

DP( 12,  8) 

1.5 

P(15,  3) 

1.2 

- 

- 

12th 

^(10,  6) 

- 

- 

- 

- 

13th 

ap(i2,  4) 

El 

- 

- 

- 

- 

Total 

44.2 

59.2 

27.9 

•Excluding  uppor-air  predictors. 
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(b)  Including  past  history 


N 

E 

0 

Predictor 

% 

rad 

Predictor 

% 

red 

Predictor 

% 

red 

A* 

mm 

AX 

29.8 

AP(  9,  5) 

7.9 

AP(10,  6) 

wBm 

AP(  8,  6) 

3.1 

P(  7,  5) 

5-3 

P(15,  7) 

1 

P(13,11) 

2.3 

AX 

3.3 

P(n#  7) 

3.4 

p(io,  2) 

1.8 

^C 

1.6 

* 

1.3 

P(12»  6) 

2.4 

P(lo,  8) 

2.0 

AP(12,  8) 

1.3 

AP(12,  4) 

1.3 

AP(17,  D 

1.4 

2P(  1,  3) 

1.1 

P(10,  5) 

1.6 

p(10,  5) 

1 .2 

- 

- 

AP(11,  5) 

1.4 

p(10,  4) 

5.1 

- 

- 

- 

- 

P(11,  5) 

1.4 

- 

- 

- 

- 

8,  4) 

1.7 

I  **3.2 

43.7 

50.9 

A* 

22.6 

AX 

25.4 

P(lo,  5) 

7.5 

P(15,  7) 

5.2 

P(13,H) 

6.3 

AP(10,  5) 

io.4 

P(»#  7) 

6.5 

*>(  8,  6) 

2.3 

P(1l,  5) 

7.4 

AX 

1.9 

P(12.  2) 

2.1 

&(  7,  3) 

2.8 

P(15,11) 

1.5 

P(15,  5) 

3.0 

- 

- 

AP(13,  9) 

1.7 

AP(12,  6) 

1.7 

- 

- 

AP(  9,  7) 

1.4 

X 

1.6 

- 

- 

4P(  7,  5) 

1.9 

- 

- 

- 

- 

♦ 

1.5 

- 

- 

- 

- 

P(l4,  4) 

1.4 

- 

- 

- 

- 

45.6 

4o.4 

27.9 
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predictors. was  the  same  as  the  first  except  the  three  past-history  predictors  were 
included.  Results  of  the  analysis  of  the  developmental  sample  for  these  two  exper¬ 
iments  are  summarised  in  Tables  6-20  and  6-21. 

The  equations  derived  in  the  two  experiments  were  applied  to  the  169  cases 
not  used  in  their  derivation.  The  results  are  presented  in  Table  6-22.  The  rms 
errors  associated  with  the  equations  using  past  history  as  predictors  are  sub¬ 
stantially  lower  than  those  not  using  past  history  for  the  two  12-hr  displacement 
predictands;  otherwise  only  slight,  if  any,  improvements  are  realised. 
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7.0  OPERATIONAL  APPLICATION  OF  THE  PREDICTION  TECHNIQUE 

Tests  on  independent  data  reveal  that  the  definition  of  best  equations  depends 
on  storm  location  and  availability  of  predictor  data.  Predictor  availability  can  be 
divided  conveniently  into  four  cases: 

Case  1.  Surface  and  upper-air  data  are  available,  but  past-history  data  are 
either  not  available  or  subject  to  large  errors  of  estimation. 

Case  2.  Surface,  upper-air,  and  past-history  data  are  all  available. 

Case  3.  Only  surface  data  are  available. 

Case  4.  Surface  and  past-history  data  are  available,  but  upper-air  data  are  not. 

In  Section  7.1,  the  recommended  winter-cyclone  equations  to  be  used  for  each  of 
the  four  data  cases  are  indicated.  A  summary  of  the  results  of  applying  each  of 
the  four  sets  of  equations  to  the  independent  data  sample  is  presented.  Section 
7.2  indicates  the  recommended  sets  of  equations  for  summer  cyclones. 

7.1  Winter  Cyclones 

Table  7-1  identifies  the  recommended  equations  for  case  1  winter  cyclones. 

All  but  two  of  the  base-technique*equations  listed  in  Section  B.  1.1.5  should  be  used: 

6  in  the  northeastern  zone  and  $  in  the  northwestern  zone.  Table  7-2  sum- 
marizes  the  results  of  applying  this  set  of  recommended  prediction  equations  to 
the  213  independent  cases. 

Thirty-one  of  the  36  base -technique  equations  (Section  B.  1.1.5)  apply  to  case 
2  winter  cyclones.  Table  7-3  identifies  the  other  five  equations  to  be  used.  These 
are  the  five  that  produced  statistically  more  significant  results  that  the  corre¬ 
sponding  base -technique  equations.  Table  7-4  is  a  summary  of  the  results  of 
applying  the  recommended  set  of  prediction  equations  to  the  213  independent  cases. 
For  case  3  winter  cyclones,  the  equations  listed  in  Section  B.1.3  should  be  used. 
k  .  For  case  4  winter  cyclones,  the  equations  listed  in  Section  B.1.4  should  be  used. 
Results  of  applying  the  recommended  case  3  and  case  4  equations  to  the 
independent  data  sample  are  summarized  in  Table  7-5. 
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Although  overall  accuracy  decreases  somewhat  as  the  number  of  available 
classes  of  predictors  decreases,  useful  predictions  appear  to  be  possible  despite 
incomplete  data. 

7.2  Summer  Cyclones 

Prediction  equations  for  case  1  summer  cyclones  are  listed  in  Section  B.2.1. 
Prediction  equations  for  case  2  summer  cyclones  are  listed  in  Section!  B.2.2. 
Prediction  equations  for  case  3  summer  cyclones  are  listed  in  8ection  B.2.3. 
Prediction  equations  for  case  4  summer  cyclones  are  listed  in  8ection  B.2.4. 
Results  of  applying  the  four  sets  of  equations  to  169  independent  cases  are 
summarized  in  Table  7-6. 
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8.0  FORECAST  EXAMPLES 


From  the  213  independent  winter  cases,  six  cyclones  have  been  chosen  to 
illustrate  the  applicability  of  certain  characteristics  of  the  forecast  equations 
to  a  variety  of  forecast  situations.  Most  of  the  forecasts  selected  were  good 
(there  is  one  poor  forecast  included);  more  representative  verification  statistics, 
based  on  the  entire  independent  sample,  can  be  found  in  Section  7.0. 

The  figures  for  this  section  are  grouped  at  the  end  of  the  section.  Each 
figure  appears  in  two  parts,  (a)  and  (b).  Part  (a)  is  the  surface  chart  at  forecast 
(initial)  time,  with  pressures  indicated  in  millibars  (hundreds'  and  thousands' 
digits  deleted).  Part  (b)  indicates  the  positions  (reading  generally  from  left  to 
right)  of  the  cyclone  12,  24,  and  36  hr  later:  circles  represent  observed  positions; 
squares  represent  positions  forecast  by  the  base  technique  on  the  unstratified  data 
sample,  using  both  surface  and  upper-air  predictors;  S' s  represent  positions  fore¬ 
cast  by  the  modified  technique  (Section  7.0),  which  uses  the  optimum  set  of  forecast 
equations;  and  triangles  represent  positions  forecast  by  equations  based  only  on 
surface  and  past-history  predictors.  Each  symbol  is  accompanied  by  the  appropriate 
central-pressure  indication. 

Tables  following  the  figures  summarize  the  pertinent  verification  statistics 
and  correspond  in  numbering  to  the  figures. 

8.1  Two  Northwestern-zone  Cyclones 

Figure  8-1  describes  a  cyclone  that,  at  initial  time  0000  GCT  18  Nov.  1959, 
was  centered  in  Saskatchewan  (54*N,  107*W)  and  had  a  central  pressure  of  1009  mb. 
The  500-mb  chart  at  the  same  time  (not  shown)  showed  a  broad  northwesterly  flow 
pattern  from  the  Canadian  Rockies  to  an  intense  trough  over  the  Great  Lakes. 

One  of  the  major  problems  presented  to  the  synoptic  forecaster  in  situations 
of  cyclones  under  a  northwesterly  flow  is  determining  if  and  when  the  cyclone  will 
"recurve."  This  cyclone  moved  rapidly  southeastward  for  the  first  24  hr,  then 
recurved  to  the  northeast.  Both  the  base  and  modified  techniques  gave  sufficient 


78 


indication  of  likely  recurvature  in  the  24-to-36-hr  interval.  The  surface-only 
forecast,  for  which  only  12-  and  24-hr  prediction  equations  were  derived,  appeared 
to  predict  recurvature  at  24  hr.  In  no  case, however,  was  the  central -pressure 
forecast  very  good. 

Figure  8-2  illustrates  the  forecast  of  a  cyclone  with  an  unusual  track.  The 
weak  cyclone  (1014  mb)  over  western  Alberta  (53*N,  117*W)  at  1200  GCT  24  Feb. 

1960  dropped  nearly  due  southward  for  36  hr,  arriving  in  central  Nevada,  850  naut 
mi  south  of  its  original  position. 

The  modified  technique  produced  the  best  forecast  for  this  case,  and,  although 
the  36-hr  displacement  error  was  fairly  large,  the  forecast  was  quite  good  qualitatively. 
The  central-pressure  forecasts  indicated  filling  whereas  in  fact  the  cyclone  deepened 
11  mb  in  36  hr. 

8.2  Two  Southeastern-zone  Cyclones 

Figure  8-3  describes  a  1002 -mb  coastal  disturbance  that  began  off  the  New 
Jersey  coast  (38*N,  74*W)  at  0000  GCT  19  Jan.  1960.  It  subsequently  moved  up  the 
coast  and  turned  eastward  in  Nova  Scotia.  At  the  end  of  the  36-hr  period,  the  cyclone 
had  deepened  24  mb,  to  978  mb. 

Both  the  base  technique  and  the  surface -only  technique  did  quite  well  with 
this  situation.  The  errors  in  displacement  and  intensification  are  small.  Even 
the  eastward  turn  was  well  handled. 

Figure  8-4(a)  is  typical  of  many  winter  weather  maps.  In  this  instance  (1200 
GCT  18  Feb.  1960),  there  is  a  deepening  cyclone  (998  mb)  on  the  Georgia- Alabama 
bolder  (38*N,  85*W).  The  500-mb  chart  (not  shown)  indicates  a  trough  over  the 
lower  Mississippi  Valley  and  a  ridge  through  New  England.  Determining  the  future 
track  of  this  class  of  cyclone  is  an  important  step  in  preparing  a  weather  forecast 
for  the  North  Atlantic  states.  This  particular  cyclone  moved  northeastward  to  the 
New  England  coast  in  36  hr,  with  moderate  deepening. 

The  displacement  forecasts  were  exceptionally  good,  with  an  error  of  less 
than  50  naut  mi  for  36  hr. 
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8.3  Two  Southwestern-zone  Cyclones 

Figure  8-5  illustrates  one  of  the  poorest  forecasts  of  the  entire  independent 
sample.  The  surface  chart  for  0000  GCT  12  Jan.  1960  shows  a  996-mb  cyclone 
over  eastern  Oregon  (44°N,  119*W).  Although  this  cyclone  moved  southeastward 
into  Arizona,  all  the  techniques  predicted  that  it  would  move  pretty  much  due 
eastward;  the  36-hr  displacement  error  was  very  large. 

Figure  8-6  represents  a  troublesome  situation  for  synoptic  forecasters. 

The  deep  cyclone  (978  mb)  over  northeastern  Colorado  (41*14,  104*W)  at  1200  GCT 
9  Feb.  1960  had  been  moving  southeastward  for  the  previous  12  hr.  The  timing 
of  recurvature  is  an  important  part  of  this  forecast.  After  12  hr  of  continued 
southeastward  movement,  the  cyclone  recurved  northeastward,  reaching  the 
Michigan-Indiana  border  after  36  hr.  The  changes  in  central  pressure  were 
small;  the  cyclone  remained  quite  intense  during  the  entire  forecast  interval. 

The  forecasts,  although  not  moving  the  cyclone  fast  enough,  did  a  good  job 
of  capturing  the  recurvature  by  indicating  a  northeastward  turn  after  12  hr. 

The  central -pressure  forecasts  were  adequate. 
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Fig.  8-2<a).  Second  cyclone-forecast  example.  Initial  conditions  at  1200  GCT  24  Feb.  1960.  Cyclone  is 
located  at  52.8*N,  116.9*W;  central  pressure  is  1014  mb. 


fig.  8-3(b).  Forecast  verification.  o  observed,  o  base  technique,  *  surface  predictors  only;  adjacent 
tends  indicate  central  pressure. 
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Fig.  8-5(b).  Forecast  verification,  o  observed,  o  base  technique,  a  surface  predictors  only,  x  modified 
technique:  adjacent  numerals  indicate  central  pressure. 
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Fig.  8 -6(b).  Forecast  verification.  o  observed,  o  base  technique,  a  surface  predictors  only,  x  modified 
technique;  adjacent  numerals  indicate  central  pressure. 


TABLE  8-1 

VERIFICATION  SUMWRY  OF  CYCLONE  SHOWN  IN  FIG.  8-1 


Forecast 

interval, 

hr 

Position-vector  error,  naut  mi 

Central-pressure  error,*  mb 

Base 

tech 

Mod  i  f  i  ed 
tech 

Surface 
data  only 

Base 

tech 

Mod i f i ed 
tech 

Surface 
data  only 

12 

12 

12 

36 

+3 

+3 

+  1 

24 

180 

156 

192 

+1 

+1 

-3 

36 

42 

60 

- 

+15 

+15 

- 

•Centra I -pressure  error  equals  predicted  minus  observed. 


TABLE  8-2 

VERIFICATION  SW*KRY  OF  CYCLONE  SHOWN  IN  FIG.  8-2 


Forecast 

interval. 

Position-vector  error,  naut  mi 

Central -pressure  error,*  mb  j 

hr 

Base 

Modified 

Surface 

Base 

Modified 

Surface 

tech 

tech 

data  only 

tech 

tech 

data  only 

12 

84 

84 

168 

+3 

+3 

0 

24 

252 

360 

432 

+9 

+9 

+9 

36 

468 

360 

- 

+19 

+19 

- 

•Centra I -pressure  error  equals  predicted  minus  observed. 


TABLE  8-3 

VERIFICATION  SUMMARY  OF  CYCLONE  SHOWN  IN  FIG.  8-3 


Forecast 
i nterva 1 , 
hr 

Posit  ion-vector  error,  naut  mi 

Central -pressure  error,*  mb 

Base 

tech 

Modified 

tech 

Surface 
data  only 

Base 

tech 

Modified 

tech 

Surface 
data  only 

12 

54 

- 

ice 

+9 

- 

+7 

24 

60 

- 

1 02 

+7 

- 

+5 

36 

158 

- 

- 

+5 

- 

- 

•Centra I -pressure  error  equals  predicted  minus  observed. 
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'  TABLE  8-4 

VERIFICATION  SUMMARY  OF  CYCLONE  SHOWN  IN  FIG.  8-4 


Forecast 

interval, 

hr 

Poslti 

on- vector  error,  naut  mi 

Central-pressure  error,*  mb  | 

Base 

tech 

Mod  1 f 1 ad 
tech 

Surface 
data  only 

Base 

tech 

Modified 

tech 

Surface 
data  only 

12 

132 

- 

48 

+10 

- 

+9 

24 

115 

- 

30 

+8 

- 

+6 

36 

24 

- 

- 

+3 

- 

■Hi 

•Centra I -pressure  error  equals  predicted  minus  observed. 


TABLE  8-5 

VERIFICATION  SUMMARY  OF  CYCLONE  SHOWN  IN  FIG.  8-5 


Forecast 

interval, 

hr 

Posit ion-vactor  error,  naut  mi 

Central -pressure  error,*  mb  j 

Basa 

Modified 

Surface 

Base 

Modified 

Surface 

tech 

tech 

data  only 

tech 

tech 

data  only 

12 

126 

144 

108 

-2 

-2 

-3 

24 

360 

360 

360 

-7 

‘7 

-10 

36 

768 

768 

- 

-8 

-8 

• 

•Centra I -pressure  error  equals  predicted  minus  observed. 


TABLE  8-6 

VERIFICATION  SUMMARY  OF  CYCLONE  SHOWN  IN  FIG.  8-6 


Forecast 

interval, 

hr 

Position-vector  error,  naut  mi 

Central -pressure  error,*  mb  j 

Base 

tech 

Modified 

tech 

Surface 
deta  only 

Base 

tech 

Modified 

tech 

Surface 
data  only 

12 

u 

136 

l8o 

-1 

-1 

+1 

24 

i9e 

228 

0 

0 

+3 

36 

19 

156 

- 

+5 

+5 

- 

•Central-pressure  error  equals  predicted  minus  observed. 
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9.0  SUMMARY  AND  CONCLUSIONS 


Equations  for  the  prediction  of  12-.  24-.  and  36-hr  displacements  and  changes 
in  central  pressure  have  been  derived  for  North  American  cyclones.  The  equations 
remained  stable  when  applied  to  an  independent  data  sample.  Comparisons  with  NWAC 
forecasts  for  the  winter  of  1962-1963  indicate  that  the  equations  yield  competitive 
results.  Equations  applicable  to  the  entire  area  considered  in  North  America  pro¬ 
duced  results  that  were  in  general  at  least  as  good  as  equations  derived  for  sub¬ 
divisions  (zones)  of  the  area. 

Equations  were  derived  from  surface  data  only,  for  application  whenever  upper- 
air  data  are  not  available.  Only  a  small  decrease  in  accuracy  was  noted  when  the 
results  of  these  equations  were  compared  with  those  derived  from  both  surface  and 
upper-air  data. 

Application  of  the  equations  is  completely  objective,  and  they  can  be  adapted 
for  use  at  a  weather  station*  or  at  a  completely  automated  weather  central. 

A  useful  by-product  of  the  data  processing  required  for  the  research  is  a 
5-yr  sample  of  error-checked  hemispheric  gridpoint  values  of  surface  pressure 
and  500-mb  heights  filed  in  a  form  suitable  for  direct  electronic  data-processing 
machine  analysis. 

The  developmental  sample  used  in  the  regression  analysis  was  derived  from 
a  series  of  computerized  objective  techniques.  This  is  in  contrast  to  the  earlier 
study  [14] ,  in  which  the  data  were  extracted  and  tabulated  manually.  Comparison 
of  results  indicates  that  the  automated  procedures,  besides  being  expeditious, 
apparently  suffered  no  information  loss  from  interpolation  errors. 

♦See  Appendix  C  for  worksheets  illustrating  manual  application. 
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10.0  RECOMMENDATIONS 

Although  the  primary  concern  in  these  studies  has  been  the  behavior  of  the 
cyclone  itself,  there  exists  a  need  to  incorporate  the  information  contained  in 
such  a  circulation  forecast  into  a  weather  forecast.  Even  a  perfect  forecast  of 
sea-level  pressure  is  of  little  practical  value  if  there  is  no  means  by  which  the 
weather  (cloudiness,  precipitation,  wind,  temperature,  etc.)  can  be  predicted  or 
inferred.  More  emphasis,  therefore,  should  be  given  to  specifying  weather  from 
circulation  forecasts  to  yield  more  useful  predictive  information. 

The  use  of  the  methods  described  in  this  study  yields  information  about  the 
sea-level  pressure  at  widely  scattered  locations.  A  more  complete  depiction  of 
the  surface  circulation  pattern  would  be  desirable.  One  possible  way  to  acquire 
a  complete  mapping  of  the  pressure  field  is  to  employ  existing  statistical  and 
dynamical  methods  to  forecast  the  entire  pressure  field  and  modify  them  with 
the  regression  predictions  of  cyclones  and  other  significant  features  such  as 
anticyclones.  Such  a  modification  could  be  accomplished  by  methods  similar  to 
those  used  in  objective  analysis  techniques. 

Thus  far,  the  moving-coordinate  experiments  have  used  only  initial-  and 
lagged-time  data.  The  incorporation  of  mid-tropospheric  circulation  prognoses 
by  numerical  weather-prediction  methods  as  possible  predictors  in  the  statistical 
analysis  appears  to  have  a  reasonable  chance  for  improving  the  results  presented 
in  this  study  and  should  be  attempted. 

The  technique  described  in  this  study  should  not  be  limited  to  extrema  in  the 
sea-level  pressure  pattern.  Other  "significant  features,"  such  as  500-mb  vorticity 
centers  or  height-change  centers,  can  be  used  as  predictands. 
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APPENDIX  A.  DATA  PROCESSING* 


One  of  the  major  problems  encountered  when  undertaking  a  statistical  problem 
of  the  nature  described  in  this  report  concerns  the  data  to  be  used  in  the  investigation. 
When  the  scope  of  the  problem  reaches  the  point  at  which  some  11,000,000  gridpoint 
values  are  extracted  manually  from  analyzed  maps,  tabulated,  and  placed  on  punched 
cards,  errors  will  persist  despite  the  most  careful  and  systematic  error-checking 
procedures.  The  best  that  can  be  hoped  for  is  a  minimization  of  ••gross"  errors. 
Although  the  data  received  for  this  project  had  already  undergone  considerable  checks 
for  errors,  it  was  felt  necessary  to  devise  additional  automatic  error-checking  pro¬ 
cedures  and  subject  the  data  to  them.  Although  only  a  small  percentage  of  additional 
errors  were  detected,  enough  were  found  to  make  the  procedures  extremely  worth 
while. 

There  were  two  basic  data  types  used  in  the  development  of  multiple  linear- 
regression  prediction  equations  for  cyclone  displacement  and  intensification  over 
North  America.  Basically,  these  are  historical  series  of  weather  maps  (predictor 
data)  and  summaries  of  cyclone  tracks  (predictand  data). 

A.l  Predictor  Data 

The  basic  predictor  data  are  a  collection  of  gridpoint  data  made  at  the 
National  Weather  Records  Center  (NWRC)  at  Asheville,  N.C.  Pressure-height  and 
temperature  at  several  levels  (including  surface,  700  mb,  and  500  mb)  over  the 
Northern  Hemisphere  were  read  from  weather  charts  analyzed  by  the  National 
Weather  Analysis  Center  (NWAC)  on  a  modified  JNWP  grid  of  1020  points.  Each 
map  was  actually  read  twice  to  eliminate  errors.  These  gridpoint  values  were 
tabulated,  placed  on  punched  cards,  and  finally  on  IBM  705  magnetic  tapes.  The 
final  set  of  data  covered  a  continuous  period  of  twice-daily  (0000-  and  1200-GCT) 
analyses  from  April  1955  through  March  1960^,  a  complete  5-yr  sample.  For  the 
three  fields  of  interest  in  significant-feature  studies  (sea-level  pressure,  700-mb 
height,  and  500-mb  height),  this  involved  more  than  11,000,000  numbers.  A  de¬ 
tailed  description  of  the  gridpoint  data-extraction  procedures  and  card-tape  for¬ 
mats  is  given  in  the  Reference  Manual  for  Climatic  Data  Computer  Tapes  (11  ]. 

A. 2  Predictand  Data 

Cyclone  tracks  were  tabulated  from  microfilmed  NWAC  surface  analyses. 

The  location  (latitude  and  longitude)  and  central  pressure  of  every  winter  (November- 
March)  cyclone  over  North  America  (as  defined  in  Fig.  2-1)  were  plotted  at 


•Most  of  the  data  preprocessing  was  cosponsored  by  the  United  States  Air  Force, 
under  Contract  AF19(626)-16,  and  the  Federal  Aviation  Agency,  under  Contract 
FAA/BRD-363. 
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APPENDIX  A.  DATA  PROCESSING* 


One  of  the  major  problems  encountered  when  undertaking  a  statistical  problem 
of  the  nature  described  in  this  report  concerns  the  data  to  be  used  in  the  investigation. 
When  the  scope  of  the  problem  reaches  the  point  at  which  some  11,000,000  gridpoint 
values  are  extracted  manually  from  analyzed  maps,  tabulated,  and  placed  on  punched 
cards,  errors  will  persist  despite  the  most  careful  and  systematic  error-checking 
procedures.  The  best  that  can  be  hoped  for  is  a  minimization  of  ••gross"  errors. 
Although  the  data  received  for  this  project  had  already  undergone  considerable  checks 
for  errors,  it  was  felt  necessary  to  devise  additional  automatic  error-checking  pro¬ 
cedures  and  subject  the  data  to  them.  Although  only  a  small  percentage  of  additional 
errors  were  detected,  enough  were  found  to  make  the  procequres  extremely  worth 
while. 


There  were  two  basic  data  types  used  in  the  development  of  multiple  linear- 
regression  prediction  equations  for  cyclone  displacement  and  intensification  over 
North  America.  Basically,  these  are  historical  series  of  weather  maps  (predictor 
data)  and  summaries  of  cyclone  tracks  (predictand  data). 

A.l  Predictor  Data 

The  basic  predictor  data  are  a  collection  of  gridpoint  data  made  at  the 
National  Weather  Records  Center  (NWRC)  at  Asheville,  N.C.  Pressure-height  and 
temperature  at  several  levels  (including  surface,  700  mb,  and  500  mb)  over  the 
Northern  Hemisphere  were  read  from  weather  charts  analyzed  by  the  National 
Weather  Analysis  Center  (NWAC)  on  a  modified  JNWP  grid  of  1020  points.  Each 
map  was  actually  read  twice  to  eliminate  errors.  These  gridpoint  values  were 
tabulated,  placed  on  punched  cards,  and  finally  on  IBM  705  magnetic  tapes.  The 
final  set  of  data  covered  a  continuous  period  of  twice-daily  (0000-  and  1200-GCT) 
analyses  from  April  1955  through  March  I960,  a  complete  5-yr  sample.  For  the 
three  fields  of  interest  in  significant-feature  studies  (sea-level  pressure,  700-mb 
height,  and  500-mb  height),  this  involved  more  than  11,000,000  numbers.  A  de¬ 
tailed  description  of  the  gridpoint  data-extraction  procedures  and  card-tape  for¬ 
mats  is  given  in  the  Reference  Manual  for  Climatic  Data  Computer  Tapes  (11  ]. 

A.2  Predictand  Data 

Cyclone  tracks  were  tabulated  from  microfilmed  NWAC  surface  analyses. 

The  location  (latitude  and  longitude)  and  central  pressure  of  every  winter  (November- 
March)  cyclone  over  North  America  (as  defined  in  Fig.  2-1)  were  plotted  at 


*Most  of  the  data  preprocessing  was  cosponsored  by  the  United  States  Air  Force, 
under  Contract  AF19(626)-16,  and  the  Federal  Aviation  Agency,  under  Contract 
FAA/BRD-363. 
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6-hr  intervals  for  the  five  winters  corresponding  to  the  gridpoint  data.  This 
included  cyclones  that  were  either  previously  or  subsequently  outside  the  area. 

Similar  procedures  were  used  for  summer  (May-September)  cyclones.  From 
these  plotted  tracks,  only  cyclones  that  were  identifiable  for  36  hr  (24  hr  for 
summer)  were  tabulated  and  put  on  punched  cards.  For  those  that  qualified, 
the  following  information  was  tabulated  for  each  synoptic  time  (0000  or  1200  GCT). 

(a)  Date  and  time. 

(b)  Latitude,  longitude,  and  central  pressure  (initial  time  tQ). 

(c)  Latitude,  longitude,  and  central  pressure  (tQ  +  12  hr). 

(d)  Latitude,  longitude,  and  central  pressure  (t^  +  24  hr). 

(e)  Latitude,  longitude,  and  central  pressure  (t^  +  36  hr). 

(f)  Latitude,  longitude,  and  central  pressure  (tQ  -  6  hr). 

The  number  of  cases  shown  in  Table  2-1  is  equivalent  to  the  number  of  cards  that 
were  punched. 

With  the  IBM  705  magnetic  tapes  and  the  punched  cards  of  cyclone  tracks 
available  for  processing,  a  variety  of  computer  programs  (Table  A-l)  was 
written  at  the  United  Aircraft  Corporation  Research  Laboratory.  These  programs 
and  their  sequential  use  are  shown  in  the  flow  diagram,  Fig.  A-l.  In  detail,  the 
programs  were  used  as  follows. 

The  first  step  was  to  process  each  of  the  Asheville  IBM  705  grid  tapes 
separately  by  the  gridpoint  data-conversion  and  -editing  program.  Since  each 
tape  contained  only  2  mo  of  grid  data,  30  tapes  were  enough  to  cover  the  full 
60-mo  (5-yr)  period.  This  program  selects  the  variables  to  be  processed: 
sea-level  pressure,  700-mb  height,  and  500-mb  height.  Other  variables  on  the 
tapes,  such  as  heights  at  other  levels  (e.g.,  300  mb)  and  temperatures  at  any 
level,  were  not  processed.  Having  made  this  selection,  the  program  then 
examined  every  "map"  (1020  gridpoints  covering  the  Northern  Hemisphere)  for 
gross  errors.  Differences  in  values  between  adjacent  gridpoints  were  computed, 
and  a  preselected  percentage  of  the  largest  differences  were  printed  out,  along 
with  appropriate  identification.  Also,  gridpoints  with  missing  data  and  values 
outside  certain  specified  limits  were  noted  and  printed  out.  This  program  also 
converted  the  selected  data  and  output  it  onto  magnetic  tape  in  a  format  acceptable 
to  the  kind  of  computer  to  be  used  for  subsequent  stages:  of  processing  (the  IBM  7090). 
The  output  of  this  program  was  a  set  of  three  7090  tapes  for  every  one  that  was 
used  as  input,  each  of  which  contained  only  information  from  one  level  (surface, 
700-mb,  or  500-mb)  for  a  2-mo  period.  Thus,  the  processing  of  the  thirty  IBM  705 
tapes  resulted  in  ninety  IBM  7090  tapes. 

The  next  step  was  very  important  and  was  strictly  a  manual  operation.  It 
consisted  of  examining  the  printout  of  missing  and  otherwise  suspicious  gridpoint 
data.  Many  of  the  data,  though  printed  out,  were  acceptable  due  to  unusually  strong 
gradients  or  intense  systems.  For  example,  all  sea-level  pressures  in  excess  of 
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Fig.  A-l.  Computer  programs  used  to  preprocess  data  for  use  in  screening  re¬ 
gression.  Numbers  appearing  within  tape  symbols  refer  to  the  number  of  magnetic 
tapes  involved. 
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1060  mb  were  printed  out.  Many  of  these  were  valid  beeause  of  the  Siberian  High, 
for  instance.  However,  there  were  some  pressures  between  1090  and  1099  mb 
that  were  actually  supposed  to  have  been  between  990  and  999  mb,  as  shown  by 
examination  of  microfilm  maps.  The  practice  then  was  to  go  from  the  printout 
to  NWAC  microfilmed  analyses  to  verify  the  validity  of  suspicious  data.  In  this 
way,  gross  errors  were  noted,  and  corrected  values  were  tabulated.  Corrections 
to  missing  data  were  also  supplied  from  examination  of  microfilmed  maps.  It 
should  be  emphasized  that  the  grid  data  received  from  NWRC  were,  on  the  whole, 
in  very  good  condition,  with  the  error  percentage  being  extremely  small;  however, 
the  nature  of  the  study  required  even  additional  scrutiny  of  the  data  to  eliminate 
so-called  gross  errors. 

The  gridpoint  corrections,  having  been  tabulated,  were  placed  on  punched 
cards  to  be  used  as  input  to  the  grid-correction  program.  All  correction  cards 
for  one  level  for  5  yr  of  one  month  were  processed  together  (e.g.,  500-mb  heights 
for  five  Januaries).  The  program  inserted  the  corrections  in  the  proper  places  and 
wrote  a  new  tape  consisting  of  the  one  level  for  5  yr  of  one  month.  Thus  the  pro¬ 
gram  not  only  inserted  the  necessary  corrections  but  also  reduced  the  number  of 
required  tapes  from  90  to  96. 

The  grid  format  used  for  the  initial  data  extraction  at  NWRC  was  a  modified 
JNWP  grid,  in  which  data  were  read  for  every  other  JNWP  gridpoint.  It  was  found 
desirable  to  have  data  at  every  JNWP  gridpoint  to  facilitate  subsequent  computations. 
The  grid-interpolation  program  accomplished  this  by  computing  values  at  the  in-be¬ 
tween  gridpoints  by  a  least-squares  fitting  of  a  quadratic  surface,  using  surrounding 
data  at  known  gridpoints.  Thus,  the  program  output  complete  maps  on  the  1977-point 
JNWP  grid.  Another  feature  of  the  program  is  that  it  used  three  tapes  as  input 
(e.g.,  five  Januaries  of  sea-level  pressure,  700-mb  height,  and  500-mb  height).  It 
read  in  the  first  record  (map)  from  each  tape,  which  were  all  for  the  same  time 
(e.g.,  0000  GCT  1  Jan.  1956),  and  output  them  as  the  first  three  records  of  the  new 
tape.  Hence,  the  final  tape  for  a  month  was  in  chronological  order  (i.e.,  three 
levels  for  0000  GCT  1  Jan.  1956,  three  levels  for  1200  GCT  1  Jan.  1956,  three  levels 
for  0000  GCT  2  Jan.  1956,  etc.,  through  31  Jan.  1960).  The  result  was  that  data 
on  the  36  tapes  using  a  modified  JNWP  grid  were  now  contained  on  12  tapes  with  a 
full  JNWP  grid,  each  tape  representing  a  complete  5-yr  month.  These  tapes  were 
then  ready  to  serve  as  the  basic  set  of  "maps"  for  predictor  information. 

The  predictand  data  were  handled  independently  of  these  programs.  The  pre- 
dictand-preprocessor  program  accepted  the  punched-card  input  of  significant-feature 
tracks  and  computed  the  two  components  of  displacement  and  the  change  in  central 
pressure,  which  were  the  predictand  forms  desired.  This  information  was  then 
placed  on  magnetic  tape.  An  important  by-product  of  this  program  is  a  printout  of 
displacements  and  changes  in  central  pressure  that  deviate  considerably  from 
climatology.  Thus,  an  error-checking  procedure  was  provided,  and  new  punched 
cards  were  prepared  when  errors  were  noted. 
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The  predictand  tape  and  the  predictor  tapes  were  then  used  as  input  to  the 
significant-featir  es  preprocessor  program.  This  program  generated  a  tape  for 
screening  regression  by  selecting  cyclones,  "reading"  maps,  and  deriving  predictors: 
the  first  cyclone  case  was  determined  from  the  cyclone-tabulation  input.  Using  the 
time  and  date  of  this  case,  the  gridpoint-data  tapes  were  searched  to  find  the  three 
data  records  corresponding  to  this  time  (one  record  for  each  map  of  sea-level 
pressure,  700-mb  height,  and  500-mb  height)  and  also  for  the  same  three  maps  12  hr 
previous  to  that  time.  These  six  records  were  read  into  the  machine.  The  next 
step  was  analogous  to  superimposing  a  moving-coordinate  grid  (centered  with  respect 
to  the  cyclone)  over  the  fixed-grid  system  and  interpolating  data  values  for  gridpoints 
in  the  moving-coordinate  system.  The  subroutine  that  computed  these  locations  is 
the  grid-transformation  routine.  Whereas  the  JNWP  or  "fixed"  grid  has  1977  grid- 
points,  the  moving  grid  has  only  221  points.  With  the  use  of  data  at  the  four  surrounding 
gridpoints  of  the  fixed  grid,  values  were  computed  by  interpolation  at  each  of  the  221 
gridpoints  of  the  moving-grid  system  for  all  six  maps  stored  in  the  machine.  The  pro¬ 
cedure  of  reading  maps  into  the  machine,  computing  locations  of  gridpoints,  and  inter¬ 
polating  data  was  accomplished  for  all  cyclones  to  be  considered.  Selection  and 
derivation  of  possible  predictors  from  this  newly  generated  set  of  data  was  then  per¬ 
formed.  Possible  predictors  consisted  of  point  values  of  sea-level  pressure,  500-mb 
height,  1000-500-mb  thickness,  and  12-hr  changes  of  these  quantities.  The  changes  in 
latitude  and  longitude  (both  in  degrees  of  latitude)  and  certral  pressure  of  the  cyclone 
were  used  as  predictands  for  12,  24,  and  36  hr. 

The  result  of  this  series  of  programs  was  magnetic  tapes  that  could  be  used 
directly  as  input  to  the  screening-regression  program  to  derive  prediction  equations 
for  cyclones  over  North  America  as  well  as  cyclones  in  Europe  and  Asia  and  anti¬ 
cyclones  in  Europe.  In  addition,  useful  by-products  include  a  relatively  error-free 
5-yr  set  of  gridpoint  data,  which  can  be  used  for  work  on  other  significant  features, 
such  as  North  American  anticyclones,  or  for  other  meteorological  problems. 
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APPENDIX  B,  PREDICTION  EQUATIONS 

The  prediction  equations  derived  from  the  regression  analysis  have  the  form 

?  A  +  AX  *■  A„X„  +  •••  ♦  A  X  ,  (13-1 ) 

0  112  2  n  n 

where  Y  is  the  predictand,  the  A's  are  constant  coefficients  derived  (by  the  method 
of  least  squares)  from  the  developmental  sample,  and  the  X's  are  the  predictors 
selected  by  the  screening  procedure. 

Each  set  of  prediction  equations  consists  of  three,  six,  or  nine  equations 
(depending  on  the  number  of  forecast  intervals):  the  three  predictands  of  northward 
displacement,  eastward  displacement,  and  change  in  central  pressure  for  the  fore¬ 
cast  intervals  of  12,  24,  and/or  36  hr. 

The  pair  of  numbers  that  is  associated  with  a  given  predictor  in  the  equations 
refers  to  the  grid  location  in  the  (k.f)-grid  system  of  Fig.  3-1.  The  symbols  used 
are  defined  in  Tables  3-1  and  3-2. 

For  operational  application,  the  equations  in  worksheet  form  represent  a 
more  convenient  means  for  solution.  An  example  of  such  worksheets  is  shown  in 
Appendix  C. 

B.l  Winter  Cyclones 

B.1.1  Equations  Using  Both  Surface  and  Upper-air  Predictors,  without  Past 
History 

B.l. 1.1  Northwestern  Zone 


N12 

-  47.1023 

+  0.0453  Z  (13,  5)  - 

■  0.031 3  Z  (9,  7)  - 

0.1212  AP  (11,7)  - 

-  0.0716  X 

-  0.0193  H 

(5,5)  -  0.0829  AP  (5,  5)  +  0.0605  P  (1 

1,3) 

£12 

-  1.1129  + 

0.0195  Z  (11,  9)  - 

0.0240  Z  (11,1)  - 

0.0321  AZ  (9,  1)  + 

0.0962  P  (11, 5) 

-  0.0876P 

(ll.D 

612 

♦  159.2940 

+  0.1743  AP  (7,  7) 

+  0.1428  AZ  (9,  5) 

-  0.5859  P  (10,  5) 

+  0.4252  P  (11,  5) 

*24 

-  11.1848 

+  0.0814  Z  (13,  5)  - 

0.0675  Z  (9,  7)  - 

0.1641 X  -  0.0464 

AZ  (5,  5) 

^24 

-  22.9553 

+  0.0414  Z  (11,  9)  - 

0.0441  Z  (11,  1)  - 

0.0728  AZ  (9,  3)  * 

0.1505  P  (11, 5) 

-  0.1244  P 

(13,1) 

624 

<  377.937 

-  0.8665  P(1 0,5)  + 

0.1752  AZ  (9,5)  ■ 

i  0.4893  P  (11, 5)  t 

0.1347  AZ  (9,  7) 
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N 


36 


12 


12 


+  17.4888  -  0.018411  (13,5)  -  0.0018  11  (7,7)  -  0.2408  X  f  0.0144  7  (1 3,  3) 

-  0.0914  Z  (9,  7)  *  0.0977  7(13,5) 

16.1877  +  0.05807.  (11, 9)  -  0.0614  7  (11,1)  -  0.0881  47(9,3)  +  0.1981  P  (11,  5) 

-  0.1783P  (13, 1) 

-  +  708.0190  -  0.5894  P  (10,  5)  +  0.1555  AZ  (9,  5)  -  0.0614  Z  (15,  1 )  +  0.1 896  AZ  (9,  7) 
B.l.1.2  Northeastern  Zone 

12.6527  «■  0.0143  Z  (15,  5)  -  0.0122  7  (7,7)  *  0.0312  7  (13,7)  -  0.0261  7  (9,7) 

*  50,1402  +  0.0327  Z  (11,9)  -  0.1557  AP  (7,  5)  »  0.0394  AZ  (1 1 ,  7)  -  0.0585  P  (15,  3) 

-  0.0148  H  (7,  7)  -  0.01437  (11,11) 

f  15 1,043  +  0-0868  A7.  (9,  3)  f  0.0606  AZ  (9,  9)  -  0.2518  P  (10,  5)  t  0.2036  P  (9,  7) 

-  0.0786  AZ  (5,  3)  -  0.0588  AZ  (15,  1 1 )  -  0.1 039  P  (15,  3) 

30.3683  +  0.0221  7  (15,5)  -  0.0461  Z  (7,  7)  ^  0.0421  7  (15,7) 

f  97.3384  +  0.0434  Z  (11, 9)  -  0.1 746  P  (1 3,  1 ) 

*  247.705  +  0.06957  (9,7)  -  0.1526  AZ  (5,  5)  +  0.1 287  AZ  (9,  9)  -  0.4580  P  (10,  5) 

-  0.083711  (13,3)  +  0.1 158  P  (15,  9)  *  0.0671  Z(7,1)  +  0.4033  AP  (1 1, 3) 

-  41.6853  +  0.0350  7(15,5)  -  0.0679  7(9,7)  +  0.0566  Z  (15,  7)  -  0.0553  AH  (9,  1) 

f  210,873  +  0.0497  Z  (13,  9)  -  0.2578  P  (13,1)  -  0.2620  AP  (7,  5)  -  0.0246H(5,9) 

*  10,99(53  f  0.1024  H  (9,  7)  -  0.1282  AZ  (5,5)  f  0.21 1 1  AZ  (9,  9)  +  0.2357  P  (15,  9) 

-  0.0854  Z  (13,  3)  -  0.2727  P  (10,  5) 

H.  1.1 .3  Southwestern  Zone 

*  163.681  «•  0.11207(13,3)  -  0.1110X  -  0.0095  Z  (7,  7)  f  0.0164  Z  (15,  7) 

0.093511  (13,3)  -  0.2717  P  (1 1, 5)  -  0.1 071  AP  (9,  7)  +  0.0418  P  (3,  9) 

+  0.2148  P  (1 1,  3)  -  0.1866  P  (1 3,  1 ) 

-  7.9962  *  0.0279  Z  (9,  9)  -  0.2068  AP  (7,  5)  -  0.024811(7,5) 

+  186.553  +  0.4135  AP  (9,  7)  -  0.1850P  (5,  5) 

-  9.6529  +  0.0298  Z  (1 3,  3)  -  0.1262X  -  0.0517  7.(9,7)  *  0.027011  (13,9) 

*  0.0289  P  (15,  3)  -  0.1503  P  (3,  3)  -  0.1432  AP  (9,  7)  *  0.0360  7.(15,5) 

-  0.1889  P(1 1,5)  ♦-  0.2569  P  (13,  3) 
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£  =•  +  18.0638  +  0.0130Z  {9, 11)  -  0.3554  AP  (7,  5)  -  0.0452H(7,5)  +  0.0636  Z  (11,  7) 

+  0.2779  AP  (13,  5)  -  0.0426  Z  (13,  3) 

f>24  =  +  202.759  +  0.4996  AP  (9,7)  -  0.0931  Z  (11, 1)  +  0.0686  H  (9,  7)  -  0.5119  P  (10,  5) 
+  0.3596  P  (11,  7)  +  0.1030  AH  (7,  3) 

N  =  +  103.849  +  0.1086  Z  (15,  3)  -  0,2766  X  -  0.0557  Z  (7,7)  -  OU747  P  (11,  9) 

36 

=  -  40.0025  +  0.0634  Z  (9,  11)  -  0.4060  AP(7,  5)  -  0.0451H(7,5) 

36 

£)_  =  +  528.061  -  0.1749  Z  (11,  1)  +  0.0770  Z  (9,  11)  -  0.3320  P(10,  5) 

36 

B.  1,1,4  Southeastern  Zone 

N  =  +  42.5547  -  0.0608  P(9,  9)  +  0.0138  Z  (13,  5)  -  0.0180  Z  (7,  7)  -  0.1205  AP (9,  5) 

1  4U 

+  0.0153  Z(13,  7) 

E  =  -  48.1421  +  0.0312  Z  (11,  9)  -  0.0356  Z  (7,  3)  -  0.1361  AP(7,  3)  +  0.0547  P(13,  7) 

1  6 

6  =  +  29.7512  +  0.4539  AP  (9,  3)  +  0.1099  AZ  (9, 5)  -  0.0467  Z  (11,  3)  +  O.OSajOH  (7,8) 

i.  Ct 

+  0.2358  AP(7,  5) 

r3  =  +  85.2001  -  0.1092  P  (9,  9)  +  0.0269  Z  (15,  5)  -  0.0212Z(7,7)  -  0.2088  AP(9, 5) 

+  0.0268  Z  (15,  7)  -  0.0185  H  (7,  9) 

£24  =  '  103-031  ‘  0-0011  Z  (13,  11)  -  0.0601  Z  (9,1)  +  0.1594  P(13,  7)  +  0.0455  Z  (11,  9) 

-  0.0374  Z  (7,  5)  +  0.0210  Z  (15,  11)  -  0.0505  AZ  (5,  3) 

0  =  +  174.012  +  0.6361  AP  (9,  3)  -  0.1564  Z  (11,  1)  +  0.0620H(7,9)  +  0.1022  AZ  (9,  9) 

24 

-  0.3513  P(5,  5)  +  0.3555  P(13,  1)  +  0.0691  AZ  (11,  7)  +  0.5116  AP(9, 5) 

-  0.3279  P(10,  5)  +  0.3313  P(ll,  7) 

5}  =  +  118.502  -  0.0443  Z  (7,  9)  +  0.0591  Z  (15,  7)  -  0.2073  AP(9, 5)  -  0.1399  P(ll,  9) 

36 

-  0.0574  AZ  (7, 5) 

£  =  -  124.854  -  0.0551  Z  (9,1)  +  0.0799  Z  (13,  9)  -  0.0779  Z  (7,  5)  +  0.2308  P(ll,  7) 

36 

-  0,0822  Z  (13,  1)  +  0.0738  Z  (3,  1) 

0  e  =  +  413.088  +  0.9645  AP(9, 3)  -  0.1444  Z(9,1)  +  0.2311  AH  (9, 9)  +  0.0670Z(9,9) 
36 

-  0.2665  P (10, 5) 
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B.  1.1.5  All  Zones 


N  =  +  83.578  3  -  0.0428X  +  0.0535  Z  (13,  5)  -  0.0286  Z  (9,  7)  -  0.0279  AP(11,  7) 

+  0.0836  AP (9,  3)  -  0.0217  H  (13,  3)  -  0.0971  P(ll,  7)  -  0.1190  AP(9,  7) 

-  0.0467  P  (5,1)  +  0.0583  P  (9,  7) 

E  =  +  6.5083  +  0.0134  Z  (11,  9)  -  0.0055  Z  (9,1)  -  0.1189  AP(7,  5)  -  0.0154  AZ  (7,  3) 

+  0.0282  AZ  (11,  7)  -  0.0226  H  (9,  3)  +  0.0314  P(13, 7)  -  0.0380  P(15,  3) 

+  0.0254  P  (5,  7)  -  0.1075  P(9, 3)  +  0.0927  P(ll,  5)  +  0.0197H(9,7) 

-  0.0117  H  (7,  7)  -  0.0157  AH  (9,  5) 

D  =  +  272.183  +  0.3942  AP(9,  5)  +  0.0590 X  -  0.5929  P(10,  5)  +  0.0550  AH  (9,  5) 

+  0.4357  P  (11,  5)  -  0.0327  H  (11,  1)  +  0.0308  Z  (9,  9)  -  0.0942  P  (5,  5) 

+  0.1863  AP(9,  7)  +  0.0459  AZ  (9,  3)  -  0.0118  Z  (13,  9) 

Nb1  =  +  20.5766  -  0.0787  X  +  0.0145  Z  (15, 5)  +  0.0108  Z  (7, 7)  -  0.0971  AP(11, 7) 

24 

+  0.0269  AZ  (13,  3)  -  0.0684  Z  (9, 7)  +  0.0263  Z  (13, 7)  -  0.0428  AZ  (7, 7) 

+  0.0362  Z  (13,5)  +  0.1026  AP  (9,  3)  -  0.0493  P  (11, 9) 

E  =  +  10.1768  +  0.0317  Z  (11,  9)  +  0.0006  Z  (9,1)  -  0.0404  AZ  (7, 3)  +  0.0914  P(13, 7) 

24 

-  0.0638  Z  (11,  3)  +  0.0348  AZ  (11,  7)  -  0.1322  AP(7,  5)  -  0.0240H(7,5) 

+  0.0326  Z  (11,  5)  +  0.0093  Z  (15,  11)  -  0.0800  P  (13,  1) 

D  =  +  706.947  +  0.108 1  Z  (11,  1)  +  0.1529X  +  0.1477  AZ  (9,  9)  +  0.1124  AP(9,  3) 

+  0.0935  H  (9, 7)  -  0.8755  P(10,  5)  +  0.0892  AZ  (9, 5)  +  0.3488  P(ll,  5) 

+  0.3652  AP  (9, 5)  -  0.0738  H  (11,  5)  -  0. 1647  P  (5, 5) 

N  =  -  9.6302  -  0.1206  X  +  0.0088  Z  (15,  5)  -  0.0676  Z  (9,  7)  +  0.0607  Z  (13, 5) 

36 

-  0.2193  P(ll,  7)  -  0.0575  AZ  (7,  7)  +  0.0619  P(15, 7)  +  0.1368  P(9,  7) 

+  0.0203  Z  (15,  7)  +  0.0369  AH  (11,5) 

E  =  -  47.2610  -  0.0091  Z  (13,  11)  -  0.0234  Z  (9, 1)  +  0.0504  Z  (11, 9)  -  0.0579  AZ  (7, 5) 
3o 

+  0.1756  P(13,  7)  -  0.0406  Z  (11,  3)  +  0.0266  Z  (15,  11)  -  0.0241H(7,5) 

-  0.1783  P(13,  1)  +  0.1276  AP  (13,  5)  r  0.0807P(3,3) 

Doc  =  +  523.605  +  0.0041  Z  (11,  1)  +  0.1900X  +  0.0596Z(9,9)  -  0.7045  P(10, 5) 

36 

+  0.1200  AZ  (9,  9)  +  0.3786  P  (13, 7)  +  0.0952  AZ  (9, 5)  -  0.  1014  Z  (13, 3) 

-  0.0909  AH  (11,  5)  -  0.1426  P(5, 7)  +  0.2990  AP(11,  7)  -  0.3200  AP(3,  3) 
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B.  1,2  Equations  Using  Both  Surface  and  Upper-air  Predictors,  with  Past 
History 


B.  1.2.1  Northwestern  Zone 

N  =  -  6.2902  +  0.7725  +  0.0121  AH  (13,  5)  +  0.0258  Z  (13,  5)  -  0.0198  Z  (9,  7) 

1 2 

-  0.0576  A 

E  -  +  0.2824  +  0.4702  AA  +  0.0169  Z  (11,  9)  -  0.0175  Z  (13,  1) 

1m 

D  =  +  173.68  +  0.0663  AZ  (9,7  )  -  0.6792  P(10,  5)  +  0. 1297  AZ  (9,  5)  +  0.5046  P(ll,  5) 
-12 

N  =  -  0.0391  +  0.6192  A0  +  0.0369  AH  (13,  5)  -  0.1532  A  +  0.0069  Z  (15,  3) 

-  0.1764  AP(13,  9)  -  0.0495  Z  (9,  7)  +  0.0495  Z  (13,  5) 

£  =  -  60.922  +  0.0343  Z  (11,  9)  -  0.0559  Z  (11,  1)  -  0.0480  AZ  (9,  11) 

24 

+  0.0953  P  (13,  7)  -  0.0460  AZ  (9,  3) 

D  =  +  347.70  +  0.1408  AZ  (9,7)  -  0.9372  P  (10,  5)  +  0. 1534  AZ  (9,  5)  +  0.5892  P(ll,  5) 
+  0.5644  APC 

N  =  -  5.4003  +  0.0638  H(13,  5)  -  0.0615  H  (7,  7) 

3b 

E  =  +  15.569  +  0.0594  Z  (11,  9)  -  0.0696  Z  (11,  1) 

ob 

=  +  597.45  -  0.477  IP  (10,  5)  +  0.2663  AZ  (9,  7)  -  0.0635  Z  (13,  1) 

3b 

B.  1,2.2  Northeastern  Zone 

N  =  -  14.451  +  O-8871A0+  0.0252  Z  (15,  7)  -  0.0168  Z  (9,7) 

X  M 

E  =  -  38.823  +  0.0213Z  (11,9)  -  0.1557  AP(7, 5)  +  0.04  05 AZ  (11,  7) 

1  M 

D  -  +  34.909  +  0.1066  AZ  (9,  3)  +  0.0705  AZ  (9,  9)  -  0.2243  P(10,  5)  +  0.1871  P(9,  7) 

1m 

N  =  -  25.394  +  1.0913  A0  +  0.055 1  Z  (15,  7)  -  0.0402  Z  (9,  7) 

24 

E  -  +  76.975  +  0.0348  Z  (11,  9)  -  0.0269  Z  (1 1,  1)  +  0.0295  Z  (13, 7)  -  0.0466  AZ  (7, 5) 

M* 

-  0.0240 H  (7,  9)  -  0.1027  P(15,  3) 

D  =  +67.160  +  0.0748  Z  (9,  7)  -  0. 1837  AZ  (5,  5)  +  0.1834  AZ  (9,  9)  -  0.1963  P (10,  5) 

24 

=  -  55.825  +  0.0400  Z  (15,  5)  -  0.07 17  Z  (9,  7)  +  0.0636  Z  (15,  7) 

36 

E  =  +  166.38  +  0,0577  Z  (13,  9)  -  0.2698  P  (13,  1) 

36 

D„.  =  +  78.475  +  0.0979 H  (9,  7)  -  0. 1847  AZ  (5,  5)  +  0.1335  AZ  (9,  9)  +  0.2659  P(15,  9) 
36 

-  0.1472  Z  (13,  3)  -  0.8691  P(10,  5)  +  0.0868  Z  (9,  3)  +  0.4643  P  (11,  7) 

+  0.5406  AP  (9, 5) 
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B.  1.2.3  Southwestern  Zone 


N  =  -  40.884  +  0.8347  A<*>  +  0.0257  Z  (15,  5)  -  0.0728 A  -  0.4314AA 

i.  ^ 

^12  =  "  7,0987  +  i-0112^*  +  0.0257  Z  (11,7)  -  0.0225  Z  (7,  5)  -  0.1399  A  P(7,  3) 

£>12  =  +  197.03  +  0.3937  AP  (9,  7)  -  0.1951  P  (5,  5) 

N24  =  +  83.301  +  0.0875  Z  (13,  3)  -  0.1685  P(ll,  7)  +  0.7803  A*  -  0.1498A  -  0.0340  Z  (9,  7) 

E  =  -52.840  +  1.5444  A  A  +  0.0275  Z  (11, 11)  -  0.0671  A  Z  (7, 3)  -  0.1696  A  P  (11, 9) 

f>24  =  +  308.01  -  1.7539  A (5  +  0.0319  Z  (9, 11)  -  0.9873  P(10, 5)  +  0.1508  A  Z  (9, 5) 

+  0.6211  P(ll,  5)  -  0.1336  AH  (11,5) 

N36  =  -  28.743  +  0.0763  Z (15, 5)  -  0.0432Z(7,9)  -  0.2483A  -  0.4141  P (11,  5) 

+  0.3510  P(13, 3)  +  0.0283  Z  (3,  9) 

Egg  =  +  256.69  +  1.5141  A  A  +  0.0142  Z  (11, 11)  -  0.0576  A  Z  (7, 3)  -  0.3116  A  P(ll,  9) 

-  0.0834  AH  (13, 3)  -  0.0739  Z  (15, 1)  +  0.1331  A  Z  (ll,  7)  +  0.0596  Z  (13,7) 

-  0.3170  P(ll,l)  +  0.43720  +  0.0272  Z  (5, 11) 

°36  =  +  361,14  +  0.0989  Z(9, 11)  -  0.0782 Z  (11, 3)  -  0.3895  P(10, 5)  -  1.6799 P(ll, 5) 

B.  1.2.4  Southeastern  Zone 

N12  =  +  51.253  +  0.5136  A 0  -  0.1049  A  Pc  -  0.0618  P(9, 9)  +  0.0189  Z  (15,  5) 

-  0.0128  H  (5,  7)  -  0.0727  A  P(ll,  7) 

^12  =  "  °'5478  +  0-5182AA  +  0.0250  Z  (11,  9)  -  0.0248Z(7,3)  -  0.0988  A  P(7,  3) 

+  0.0304  A  Z  (11,  7) 

=  +  58.190  +  0.4228A  P(9,3)  +  0.0970  A  Z  (9, 5)  -  0.05dlZ(U,3)  +  0.0317H(7,9) 

+  0.3154  A  P(7,  5)  +  0.0471  AH  (9,9)  +  0.0782  P(13,  5)  -  0.1614  P(5, 5) 

-  0.8648  Ac  -  0.4314  P(10,  5)  +  0.5014  P(ll,  5) 

N24  =  -  2.0672  +  1.2065  Ac  -  0.0375  Z  (9,  9)  +  0.0399 H  (13,  7)  -  0.0394  A  Z  (7, 5) 

E24  =  -  109.23  +  0.7675  A  A  +  0.0156  Z  (13,  11)  -  0.0595  Z  (7,  3)  +  0.1301  P  (13,  7) 

+  0.0317  Z(1 1,9)  -  0.0526  Z  (15,  1)  +  0.0544  A  Z  (11,  7)  +  0.0520  Z  (3,1) 

D24  =  +  313.43  +  0.6608  A  P(9,  3)  -  0.1616  Z  (11,  1)  +  0.0555H(7,9)  +  0.1216  A  H  (9,  9) 

-  0.3220  P(5,  5)  +  0.3255  Ap(ll,  9)  +  0.1881  P(13,l)  ♦  0.5986  A  P(9, 5) 

-  0.8641  P(10,  5)  +  0.8879  P(ll,  5) 
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N36  =  +115.72  -  0.0413  Z  (7,  9)  +  0.0571  Z  (15,  7)  -  0.2036  A  P  (10,  5)  -  0.0708  A  Z  (7,  5) 

-  0.1397  P(ll,  9) 

E36  =  -  85.596  +  1.1994  A\  -  0.0510  Z  (9,1)  +  0.0658  Z  (13, .9)  -  0.0739  Z  (7,  5) 

-  0.0799  Z(  15,  1)  +  0.1902  P(1 1,7)  +  0.0794  Z  (3,  1) 

D36  =  +  393.97  +  0.7571  A  P(9,  3)  -  0.04  7  0  Z  (9,1)  +  0.2290 A  H  (9,  9)  +  0.0689  Z  (9,  9) 

-  0.6165  P(10,  5)  -  0.0729H  (13,5)  -  0.1287  A  Z  (3,  7)  +  0.5942A  P(9, 5) 

+  0.3111  P(13,  7) 

B.  1.2.5  All  Zones 

Nj2  =  -  68.994  +  0.7884  A0  +  0.0138  Z  (15,  5)  -  0.0266  Z  (9, 7)  +  0.0197  H  (11,  7) 

+  0.0582  P (13, 3)  -  0.0250 X 

Ej2  =  -  26.016  +  0.5644  AX  +  0.0184  Z  (11,  9)  -  0.0061  Z  (9,1)  -  0.0132A  Z  (7,  3) 

+  0.0260 A  Z  (11,  7)  +  0.0870  P(13,  5)  -  0.0582  P(15,  3)  -  0.0144  Z  (9, 3) 

-  0.0169 A  Z  (7, 5) 

Di2  =  +  101.86  +  0.3998  A  P (9, 5)  +  0.0578  X  -  0.0295  Z  (11,  3)  +  0.0422  H  (9, 5) 

+  0.0514  A  Z  (9, 9)  +  0.088 1 A  Z  (9, 3)  -  0.6841  P(10,  5)  +  0.6152  P(ll,  5) 

-  0.6082  A0  -  0.0340 H  (7,  1) 

N24  =  -  108.07  +  0.8390  A*  +  0.0363  Z  (15, 5)  -  0.0574  Z  (9, 7)  +  0.0339H  (11,  7) 

-  0.0481  X  +  0.0884  P(13,  3)  -  0.0253  A  Z  (7,  7) 

£24  =  -  125.19  +  0.7555  AX  +  0.0087  Z  (13,  11)  -  0.0074  Z  (9,1)  +  0.0207  Z  (11,  9) 

-  0.0553  A  Z  (7,  3)  +  0.1308  P(13,  7)  -  0.0451  Z  (11,  3)  +  0.0479  A  Z  ( 11, 7) 

♦  0.0185  H  (13,  7)  -  0.0911  A  P(9,  9) 

f>24  =  +  324.11  +  0.141SU  -  0.0658  Z  (11,  3)  +  0.0599  Z  (9, 7)  +  0.1114  A  Z  (9,  9) 

-  1.3105  P(10,  5)  +  0.5747  A  Pfi  +  0.8568  P(  11, 5)  +  0.3902 A  P (10, 5) 

-  0.2232  P  (5, 5)  +  0.3506  P  (9,  3)  -  0.8 131  A$ 

ft36  =  -  114.17  +  0.817&A*  +  0.0361  Z  (15,  3)  -  0.0300  Z  (9, 9)  +  0.0503  H  (13,  7) 

-  0.0869 X  -  0.0388  Z  (9, 7)  +  0.0891  P(15,  7)  -  0.0426  AH  (7,  7) 
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E  .  =  -  27.865  +  1.0232  AA  -  0.0043  Z(13,  11)  -  0.0359  Z(9,  1)  +  0.1710  P(13,  7) 

36 

-  0.2020  P(13,  1)  -  0.0256  H(7,  5)  +  0.0195  Z{5,  1)  +  0.0319  Z(ll,  9)  -  0.0242  AZ(7,  3) 

+  0.0231  Z(15, 11)  -  0.0310  Z(15,  1)  -  0.0410A  +  0.0984P(3,3)  +  0.0574  AZ(11,  7) 

-  0.1694  AP(9,  9)  +  0.1183  P(7, 9)  -  0.1189  P(9,  7) 

=  +  519.82  -  0.0530  Z(ll,  1)  +  0.0251  Z(9,  9)  +  0.1879  A  -  1.0954  P(10,  5) 

36 

+  0.4614  P(ll,  7)  +  0.6281  AP(10,  5)  -  0.0824  H(13,  3)  -  0.1498P(5,  7) 

+  0.1414  AZ<9, 9)  +  0.0705  H  (9,  7)  +  0.3240  P(9,  3) 

B.1.3  Equations  Using  Surface  Predictors  Only,  without  Past  History  (All  Zones)* 

N  =  +  31.510  -  0.0646  A  +  0.0686  P(15,  7)  -  0.0197  P(10,  6)  -  0.1428  AP(9,  7) 

1  U 

+  0.0881  AP(7,  3)  -  0.1153P(11, 7)  -  0.0635 P(6,  2)  +  0.1683 P(13,  3)  +  0.0938  AP(10,  2) 

+  0.0752  AP(11, 9)  +  0.0209  P(3,  9)  -  0.0839  P(13, 1) 

Ej2  =  -  0.9259  +  0.0297  P(13, 11)  -  0.0921  AP(7,  5)  -  0.0327  P(ll,  1)  +  0.0829  P(6,  6) 

+  0.0106 P(13,  7)  -  0.0347  P(17,  3)  -  0.2225 P(9,  3)  +  0.1426  AP(10, 4)  +  0.0786P(7, 1) 

+  0.1944  P(13,  5)  -  0.1091  P(14,  4)  -  0.0729  AP(8,  6)  -  0.0614  AP  (3,  5) 

Dj2  =  +  185.54  +  0.0587  AP(9,  5)  +  0.0845A  -  0.6215  P(10,  5)  +  0.5750  P(ll,  5) 

-  0.1925  AP(5  ,  3)  +  0.3847  AP(10,  4)  +  0.1764  AP(10,  8)  -  0.6411  P(17,  5) 

-  0.0840  P(6,  6)  +  0.1846  AP(8,  6) 

N24  =  +  34.478  -  0.1274  A  +  0.0424  P(17,  7)  +  0.0807  P(10,  6)  +  0.1222  P  (14,  4) 

-  0.1517  P(7, 1)  +  0.1925  AP(7, 1)  -  0.2163  AP(9,  7)  -  0.2590  P(ll,  7)  +  0.1480  AP(11, 9) 

+  0.0534 P(l,  7)  +  0.0912 P(15,  7) 

E  =  +  92.505  +  0.0678  P(13, 11)  -  0.1242  P(13, 1)  +  0.1995  P(13., 7)  -  0.0674  A  +  0.1562  <f> 

-  0.1682  AP(7,  5)  +  0.1053  P(5,  7)  -  0.1779  P(10,  2)  +  0.1704  AP(11,  5) 

-  0.0747  P(17,  3)  -  0.0942  P(9,  7) 

D  =  +  295.45  +  0.2037  A  +  0.0601AP(9,  3)  -  0.4234  AP(5,  3)  -  0.8318  P(10,  5) 

+  0.4815  AP(10,  6)  +  0.4700  P(12,  6)  -  0.2154  P(17,  3)  +  0.5081  AP(10,  4)  +  0.1,588  P(ll,  1) 

B.1.4  Equations  Using  Surface  Predictors  Only,  with  Past  History  (All  Zoneji)* 

N  -  -  4.3017  +  0.9612  A0  -  0.0359  A  +  0.0148  P(17,  7)  -  0.1004  P(ll,  7) 

X  c* 

+  0.0458  P(14  ,6)  -  0.0572  <P  +  0.0499  P(13 ,3)  +  0.0520  AP(9  ,3) 

♦Twelve-  and  24-hr  predictions  only. 
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E12  =  +13.020  +  0.6839A \  +  0.0263P  (13, 11)  -  0.1078  P  (10,  2)  +  0.1472  P  (12,  6) 

-  0.0201  X  -  0.0343  P  (15,  3)  -  0.0823  A  P  (8,  6)  +  0.0485*  +  0.0366P(5,7) 

+  0.0887  A  P  (11, 5)  -  0.0834  P  (10,  6) 

D12  =  +  87.703  +  0.1519  A  P  (9,  5)  +  0.0707  X  -  0.7367  P  (10,  5)  +  0.6389  P  (11,  5) 

-  0A822AP(5,  3)  +  0.3093  A  P  (10, 4)  -  0.8537  A*  +  0.2994  AP£ 

N  =  +  7.7605  +  1.2399  A*  -  0.0683X  +  0.0813 P  (17,  7)  -  0.1776  P  (11,  7) 

+  0.0646AP(15,  5)  -  0.1017«#»  +  (Lr002P  (14, 4)  +  0.1133 AP  (11,  9) 

-  0.1203A  P(10,  6)  +  0.1161  AP(11,  3) 

E24  =  -  0.6236  +  1.1414  AX  +  0.0678  P  (13, 11)  -  0.1650P(13, 1)  -  0.0496X  +  0.1121* 

+  0.1400  P  (13,  7)  +  0.0792  P  (3,  3)  -  0.0564  P  (17,  3)  -  0.1295P(9,  3) 

+  0.1230  A  P  (13,  5)  -  0.1124  AP(8,  6)  +  0.0585P(6,8) 

=  +  368.81  +  0.1665  X  +  0.2518  A  P  (9  ,  3)  -  0.3945  AP  (5,  3)  -  1.1104  P (10,  5) 

+  0.5362  A  P^  +  0.6422  P  (11, 5)  -»  0.4337  A  P  (10,  5)  -  0.99  34A*  +  0.2304  P  (11,  7) 

-  0.1501  P  (17,  3) 

B.1.5  Equations  Using  Derived  Predictors,  without  Past  HiBtory  (All  Zones)* 

N24  =  +  77.026  +  0.5702^  -  0.0725X  +  0.0388  P(10,  ,6) >  O.016*AH 6) 

+  0.6157  AU?  -  0.1218  P  (10,  6)  +  0.1569£  T  (8,  8)  +  0.0499  Z  (14,  4) 

-  0.0166 H  (6,  8)  -  0.0493H?(12,  2) 

t  =  +  8.2071  -  0.6142  u  0.0205  Z  (12,  10)  +  0.0158A  (9,7)  -  0.042BA  Z  (6,  5) 

7  H5 

-  0.0248Z  (11,  1)  +  0.4823 -  0. 1760  A  ^ (10, 9)  +  0.0610 11  (5,  5) 

-  0.0640  P  (17,  3)  +  0.0355  AfT  (13,4) 

N  =  +  457.88  -  0.5449 fV  |  +  0.1336X  +  Q.1093A  Z  (9, 8)  -  0.4647 P  (10,  5) 

+  0.0758  A  Z  (10,  3)  +  0.0356A„  (6,4)  -  0.1696  A_,  (9, 10)  +  0.0513A.,  (9,7) 

H  5  r  rf 

+  0.4407  A  P  (9,  6)  -  0.8333  v'  +  0.2719  AP  (10,  3) 

B.1,6  Equations  Using  Derived  Predictors,  with  Past  History  (All  Zones)* 

N  =  +  100.07  +  0.9163  A*  +  0.0559  Z  (14,  4)  -  0.0119Z(7,7)  -  0.0532X  +  0.1685£  (8,8) 

x 

-  0.1192 P (10,  6)  -  0.0588 H? (12,  2)  +  0.4108  Au?  +  0.1687  I  V  f  +  0.1182^(14,  2) 
•Twenty-four-hour  predictions  only. 
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E24  =  +  101.34  +  0.9242  AX  -  0.2744U  +  0.0157 Z  (12,  10)  +  0.1687v?  -  0.0045H  (6,8) 

-  0. 1505AE  (10, 9)  -  0.1523 P  (13, 1)  +  0.0l29Au  (9,  7)  -  0.0369  H_  (12,  2) 

l  H  7 

+  0.7737  v  +  0.0640  P  (5,  9)  -  0.4488  v'  +  0.0196  A„„  (6,  9) 

7  7  H7 

D24  =  +  459.74  -  0.5237  |  V!  +  0.1378  X  +  0.1146  A  Z  (9,  8)  -  0.4671  P  (10,  5) 

+  0.0730AZ  (10  ,  3)  +  0.0348  A  (6,4)  +  0.1851  AP  -  0.1660  A  (9,10) 

H  C  b 

+  0.0490 A  (9,7)  +  0.3887  AP(9,  6)  -  0.7968v^_  +  0.2479  A P (10,  3) 

H  7 

B.2  Summer  Cyclones 

B.2.1  Equations  Using  Both  Surface  and  Upper-air  Predictors,  without  Past 
History  (All  Zones)* 

N12  =  +  143.86  -  0.1334  AP  (9,  7)  -  0.0298  P  (15,  7)  -  0.0579  P(ll,7)  +  0.0093  AP(7,  5) 

+  0.0025  AP  (11,  9)  -  Q.0194!X  -  0.0442  P  (9,1)  -  0.0197  Z  (9,  7)  4  0.0550  Z  (13,  5) 

-  0.0207H(13,  3)  +  0.0424  P (15,  9)  -  0.0157H(9,5)  -  0.0500P(5, 1) 

-  0.1096  AP  (10,  5)  +  0.0875  AP  (11,  3) 

E12  =  +  64.417  4  0.0002  Z  (11, 9)  -  0.0358  Z  (9,  3)  4  0.0448  Z  (11,  7)  -  0.0201  H  (13,  5) 

4  0.0074  Z  (3,  7)  -  0.0588  P  (11, 1)  4  0.0724  AP  (11,  7)  -  0.0757  AP  (9,  5) 

4  0.0240  AZ  (11,  3)  -  0.0164  AZ  (7,  5) 

°12  =  +  68  976  +  0.2442  AP  (9,  5)  -  0.5834  P  (10,  5)  4  0.5122  P  (11,  5)  -  0.1222  AP  (7,  5) 

4  0.2946  AP (10,  5)  -  0.1762  AP(7,  3) 

N24  =  4  223.78  4  0.0497  P(15,  7)  -  0.1484  P(ll,  7)  -  0.2225  AP (9,  7)  4  0.1030  AP (7,  5) 

-  0.0363X  -  0.1190  P  (9,1)  -  0.0501  H  (9,  7)  4  0.0817  Z  (13,  5)  4  0.0147  Z  (15, 11) 

-  0.0478  Z  (7,1)  4  0.0306 H  (5, 1)  -  0.1760  AP (10,  5)  -  0.0292H  (13,  3) 

+  0.1469 AP(11,  3) 

E24  =  +  106.68  4  0.0181  Z(ll, 9)  -  0.0810Z(9,3)  4  0.0475  Z  (11,  7)  4  0.0193Z(3,7) 

-  0.0342H(13,  5)  -  0.1023P(15,  3)  4  0.0273 Z  (11,  5)  4  0.0489  AZ (11,  7) 

-  0.0359  AZ  (7,  5) 

D  =  4  189.8  -  0.9737  P  (10,  5)  +  0.6694  AP  (10,  5)  4  0.7810  P  (11,  5)  -  0.2985  AP  (7,  3) 
•Twelve-  and  24-hr  predictions  only. 
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B.2.2  Equations  Using  Both  Surface  and  Upper-air  Predictors,  with  Past  History 
(All  Zones)* 

N  t  79.L51  +  0.41 10  A0  *  0.0165  P  (15,  7)  -  0.0864  P  (11 , 7)  -  0.1686AP  (9,  7) 

1  Ci 

-  0. 1904  AX  +  0.0436  AP  (1 3,  9)  -  0.0166X  -  0.0790P(9,1)  +  0.0717  P  (1 3,  3) 

t-  0.0601  AP(1 1,3)  -  0.0228  H  (9,  7)  t  0.0226  Z  (13,  7) 

E  -  t  8.8369  +  0.4826  AX  t  0.0034  Z  (1 1 , 9)  -  0.0299  Z  (9,  3)  +  0.0321  Z  (11,  7) 

1  & 

-  0.0649  AP  (9,  5)  +  0.0802  AP  (13,  5)  -  0.0656  AP  (9,  3)  +  0.0098  AZ  (3,  7) 

-  0.010811(13,5) 

D  =  +  55.373  *  0. 2200  AP  (9,  5)  +  0.4158AX  -  0.5600  P  (10,  5)  +  0.5028  P  (11 , 5) 

1  A 

-  0.2317  AP(7,  3)  +  0.2558  AP  (10,  5) 

N  =  ♦  106.00  f  0.4652  A0  ♦  0.0565  P  (15,  7)  -  0.0594  P  (11 , 7)  -  0.0497  Z  (9,  7) 

+  0.0601  Z(13,  5)  +  0.0146Z  (15,  11)  -  0.08690  -  0.0991  P(5,1)  -  0.0247H(9,5) 

-  0.131 1AP  (11, 7) 

E24  -  +  4.5897  *  0.5862  AX  +  0.021 5  Z  (1 1 ,  9)  -  0.0623Z(9,  3)  +  0.0375  Z  (11,  7) 

-  0.1754  AP(9,  3)  +  0.1 143  AP  (13,  5)  +  0.0457  AZ  (11,  7)  -  0.0326  AZ  (7,  5) 

D  =  +  216.82  -  0.9411  P(10,  5)  +  0.4080  AP  (10,  5)  +  0.6747  P  (11, 5)  -  0.1632  AP  (7,  3) 

+  0.0688  P  (l  1, 9)  +  0.1 100  AP  (15,  7)  -  O.2144A0  +  0.0150H(7,7)  -  0.0262  H  (11,  1) 

+  0. 1169  AZ  (9,  5)  -  0.1010  AZ  (5,  1 ) 

B.2,3  Equations  Using  Surface  Predictors  Only,  without  Past  History  (All  Zones)* 

N  =  +  33.647  -  0. 1 288  AP  (10,  G)  +  0.0830  P  (15,  7)  -  0.1 150  P(ll ,  7)  +  0.1165  AP  (12,  8) 

X  u 

-  0.1624  AP  (9,  7)  -  0.0288  X  -  0.1226P(9,1)  +  0.0918  AP  (8, 4)  +  0.0896  P  (11,  3) 

+  0.0346  P (15, 11) 

E  =  +  48.628  +  0.0161  P(12,  8)  -  0.0895  P (11, 1)  -  0.1745  AP  (8,  6)  +  0.1988  AP (11, 5) 

X  £t 

-  0.1271 P (10, 5)  +  0.1493  P (12,  6) 

D  =  +  92.725  +  0.5253  AP  (9,  5)  -  0.1191  P(7,  5)  +  0.0518 P (10,  8)  -  0.5604 P (10,  5) 

X  m 

+  0.5332  P  (11, 5)  -  0.29  8  0  AP  (8,  4)  -  0.1978  AP  (15,  3)  +  0.1694  AP  (10,  4) 

N  =  -  0.6180  +  0.1330 P ( 15, 7 )  -  0.1753P(10,  6)  -  0.0611  X  -  0.2532 AP (9,  7) 

+  0.1694  AP(7,  5)  +  0.0763  P  (15, 11)  -*•  0.2037  P  (14,  4)  -  0.2459  P(12,  6)  -  0.2e64P(9,  1) 
+  0.2411  P  (9,  5)  ^  0.1613  AP  (12,  8)  -  0.2495  AP  (10,  6)  +  0.1659  AP  (12,  4) 

♦Twelve-  and  24-hr  predictions  only. 
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E  =  -  40.365  +  0.1118P(13,  9)  -  0.2302  P  (12,  2)  -  0.3155  AP  (8,  6)  +  0.1695  AP  (11,  7) 

+  0.3252  P  (13,  5)  -  0.0635  X  +  0.0620  P  (1,7)  +  0.1064  P  (5,  3)  +  0.2779  AP  (11,  5) 

-  0.2090 P (10, 4)  -  0.1263 P (15,  3) 

D  =  +  189.80  -  0.9737  P  (10,  5)  +  0.6694  AP  (10,  5)  +  0.7810P  (11, 5)  -  0.2985  AP(7,  3) 

B.2.4  Equations  Using  Surface  Predictors  Only,  with  Past  History  (All  Zones)* 

N  =  +  26.565  +  0.6489  A0  -  0.1402  AP  (10,  6)  +  0.0700  P  (15,  7)  -  0.0942  P  (11,  7) 

1  Li 

-  0.0434  4>  +  0.0771  AP  (12,  8  )  -  0.0982  AP  (1,  3) 

E  =  +  8.3190  +  0.6092  AX  -  0. 1325  AP  (8,  6)  +  0.0222  P  (13, 11)  -  0.0586P  (10,  2) 

X  Li 

+  0.1179  P  (12,  6)  +  0.0140  AP  (12,  4)  -  0.0922  P  (10,  5)  +  0.1231  AP  (11,  5) 

D  =  +  83.916  +  0.6191  AP(9,  5)  -  0.1870P(7,5)  +  0.4597AX  -  0.1184AP  +  0.0626  P  (10,  8) 

X  Li  C 

-  0.2804  AP (17, 1)  -  0.5531  P (10,  5)  +  0.2779  P (10,  4)  +  0.3165P(11,  5) 

-  0.2237  AP  (8, 4) 

Nn>1  =  -  75.979  +  0.8567  A0  +  0.0990  P  (15,  7)  -  0.2140  P  (11,  7)  -  0.2568  AX 
24 

+  0.0824  P  (15, 11)  +  0.1295  AP  (13, 9)  -  0.2187  AP  (9,  7)  +  0.1684  AP  (7,  5) 

-  0.0799  0  +  0.1113  P  (14, 4) 

E„  .  =  +  11.630  +  1.0631  AX  +  0.1002P(13, 11)  -  0.1772  AP(8,  6)  -  0.2987 P (12,  2) 

24 

+  0.1867  P  (13,  5)  +  0.1556  AP (12,  6)  -  0.0395X 
D  =  +  189.80  -  0.9737  P  (10,  5)  +  0.6694  AP  (10,  5)  +  0.7810  P  (11, 5)  -  0.2985  AP  (7,  3) 


♦Twelve-  and  24-hr  predictions  only. 
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APPENDIX  C.  PROCEDURE  FOR  MANUAL  OPERATIONAL  PREDICTION 


C.l  Construction  of.  Worksheets 

The  following  example  demonstrates  the  procedure  for  making  an  operational 
forecast  of  the  displacement  and  change  in  central  pressure  of  winter  cyclones  for  12, 

24,  and  36  hr.  The  case  to  be  illustrated  is  hypothetical,  but  gives  a  general  idea  of 
how  worksheets  for  other  equations  can  be  derived.  Although  other  equations  may  be 
recommended  for  specific  situations  (Section  7.0),  the  unstratified  (base-technique)  set 
of  equations  (Section  B. 1.1.5)  was  chosen  for  this  example. 

Besides  arranging  the  predictors  on  the  worksheets  in  a  form  convenient  for 
tabulation  and  computation,  the  following  modifications  are  desirable  in  preparing  work¬ 
sheets. 

(a)  Transform  the  predictands  so  that  the  final  result  of  the  computation  yields 
the  forecast  latitude,  longitude,  and  central  pressure.  ^This  can  be  done  easily  by 
adding  the  initial  latitude  and  central  pressure  to  the  N-  and  D-equations  of  Appendix  B, 
which  in  effect  is  making  these  initial  conditions  predictors  multiplied  by  a  coefficient 
of  1.  If  an  initial  condition  is  already  a  predictor,  its  coefficient  is  modified  by  adding 
1.000  to  it.  Thus  fort)  (12-hr  prediction  of  central  pressure),  P(10,5),  which  has  a 
coefficient  of  -0.5929  in  Section  B.  1.1.5,  becomes  0.4071.  The  longitude  transformation 
is  a  little  more  involved.  There,  it  is  necessary  to  convert  the  result  from  degrees  of 
latitude  to  degrees  of  longitude  before  adding  on  the  initial  longitude.  Thus  the  final 
longitude  can  be  expressed  by 

<t>  +  <PF 

Forecast  longitude  =  \  +  E  sec  — - - ,  (C-l) 

A 

where  \  is  the  initial  longitude  in  degrees  of  longitude  (negative  if  east  longitude),  E 

is  the  predicted  eastward  displacement  in  degrees  of  latitude,  <t>  is  the  initial  latitude  in 

degrees  of  latitude,  and  <j>  is  the  forecast  latitude  in  degrees  of  latitude. 

r 

(b)  Transform  the  500-mb  height  and  1000-500-mb  thickness  predictors  and 
their  time  changes  to  the  recently  adopted  decameter  system  by  multiplying  each  by 
3.2808. 

(c)  When  a  1000-500-mb  thickness  chart  is  unavailable,  the  predictor  can  be 
computed  accurately  enough  by  concerting  the  sea-level  pressure  to  1000-mb  height  by 
use  of  the  hypsometric  formula.  If  a  mean  temperature  of  0*C  between  the  surface  and 
1000  mb  is  assumed  (which  is  reasonable  for  winter  situations),  the  following  equation 
can  be  used: 


H  =  Z  -  0.8  (P  -  1000),  (C-2) 

where  H  is  the  1000-500-mb  thickness  in  decameters,  Z  is  the  500-mb  height  in 
decameters,  and  P  is  the  sea- level  pressure  in  millibars.  Base-  technique  worksheets, 
shown  in  Figs.  C-l  through  C-4,aref<r  UBe  mIkd  pessure  data  are  in  millibars  and  height 
data  are  in  decameters  (Dm). 
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NORTH  AMERICAN  CTCLOtC  PREDICTORS 
(RASE  TECHNIQUE) 


‘o' 

SiifftM  Prtdirtwi  iMd 

t  AF 


Ut(0)  M2. 3 

Lon*U)  go,  2 
P(1 1 ,9)  1014 


196_ 


iSSsJL  Pradlctora  (1*) 
t  AZ 


1 )  £21 
2(15,11)  £2  £ 


2(9,9)  £2  3 

P(5,7)  iftl£  2(11,9)  £27 

P(7,7)  1017  2(13,9)  £29 

P(9,7)  10)3  1013  Q  2(7,7)  £39 

P(1I,7)  1013  )Ol£  -2  2(9,7)  £3£ 

**(  1 3»7 )  1QP4  2(ti ,7)  £3  8 

P<  1 5,7 )  1015  2(13,7)  £M2 

P<5,5)  1018  2(15,7)  £M£ 

P(7,5)  1015  AO  11  M  *.(7,5)  £98 

P(9,5)  1007  1009  -2  *.(9,5)  £H3 

P(tO,5)  1001  1(11,5)  £99 

P(n,5)  1QQ7  1019  ~7  *i(i3,5)  £££ 

*(3.3)  ifiJJL  -1013  . O....  2(15,5)  ££3 

P(9,3)  toil  1012  -1  2(7,3)  £41 

P(13.3)  ion  Z(9,3)  ££9 

P<  *  5,3 )  1019  Z(i  1 ,3)  £43 

P(5,l )  1017  *1(13,3)  £47 

P(11,1)  IQ  IS  2(  9,1)  £72 

P('3,1 )  <017  2(11,1)  £79 

Thieknaaa  (1000-500  b^)  Pradlctora 

H(7,7)  -  2(7,7)  -  0.8  ^>(7,7)  -  1000  7  -  £2£ 

H(9,7)  -  2(9,7)  -  0.8  />(9,7)  -  1000_7  -  £2M 
H(7,5)  -  2(7,5)  -  0.8  />( 7,5)  -  1000_7  -  £3  4 
H(11,5)  -  2(11,5)  -  0.8  if>(l1,5)  -  1000  7-  £93 
8(9,3)  -  2(9,3)  -  0.8  ^>(9,3)  -  1000_7  -  ££Q 
Hd3,3)  -  2(13,3)  -  0.8  Z"P(13.3)  -  1000_7  -  ££H 
M(11,1 )  -  2(11,1)  -  0.8  £>(11,1)  -  1000_7  -  £42 

*8(9,5)  -  A2(9,5)  -  0.8  Z"*P(9,5 )J  -  -W.  4 
*8(11,5)  •  *2(11,5)  .  0.8  r*P(l1.5)J7-  1.7 

Fig,  C-l.  Data-tabulation  worksheet. 


£2  8 

JL3L1 
£3  9 


£99 

££Q 

££3 


*54*4 

£43 

■*£4.6 
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NORTH  AMERICAN  CYCLONE  PREDICTION  (12  HOUR) 


lxtituue 

ut(0)  W1.3  X  1.0000 

LongU)  30  2  X  -0.0t28 

P(9,7)  /OI3  X  0.0583 

P( 1 1 ,7)  10  13  X  .0.0971 

P(5,1 )  /Oil  X  .0.0.67 

AP(9,7)  O  X  .0.1190 

AP(11,7)  -2  X  .0.0279 

AP(9,3)  -  1  X  0.0836 

4(9,7)  3^3  3  X  .0.0938 

4(13,5)  SSS  X  0.1755 

H03.3)  SSH  X  .0.0712 

SUM  OF  PRODUCTS  -w-fcll 
COUSTaMT  ADDITIVE  83.5783 
FORECAST  LATITUDE  W3.Q 

CENTRAL  PRESSURE 


longitude 


P(5,7) 

IQ  UL 

X  0.0254 

P(  1 3  »7 ) 

)Pl.fe 

x  0.0314 

P(1 1 ,5) 

10. CJ. 

X  0.0927 

P(9,3) 

'Oil 

X  .0.1075 

P(15,3) 

JOii- 

X  -O.O38O 

6P(?,5) 

H 

X  .0.1189 

4(11,9) 

JlAJ— 

x  0.0440 

4(9,1) 

-S.7.A-, 

X  -0.0180 

64(11,7) 

_=J _ 

X  0.0925 

44(7,3) 

_ 

x  .0.0505 

H(7,7) 

■53.5 

X  .O.C384 

H(9,7) 

,£Ad.. 

X  0.0646 

H(9,3) 

X  -0.0741 

6H(9,5) 

2±Li. 

x  -0.0515 

SUM  OF 

PRODUCTS  ~IO  ~lS 

CONSTANT  ADDITIVE  6.5083 


Long(l)  <10.  JL  X  0.0590 

P(5,5)  tO  1  g  X  -0.C942 

P(lO,5)  IOOI  X  0.*071 

P(  1 1 ,5)  1007  X  0.4357 

AP(9,7)  Q  X  O.I863 

6P(9,5)  -3.  X  0.3942 

41(9.9)  -S'.!*  X  0.1010 

4(13,9)  513  X  -0.0387 

*4(9,3)  ~  H  X  0.1506 

H(ll.t)  All  X  -0.1073 

6H'9,5)  -S.fe  X  0.1804 

SUM  OF  PROOUCTS  JJJt&A 
CONSTANT  ADDITIVE  272.183 
FCST  CENTRAL  PRESSURE  m, 


(1)  Toui-  -H.aHJ _ 

(2)  1/2  (initial  ♦  feet  let.)  H2.7 

(1)  Sec  (2)  13  47  X  (1)  -  ~  5.8 

(4)  Initial  Long.  U)  _ *}Q.  3. 

(5)  FORECAST  LONG.  £{})  *  U)  7VH.H 


FORECAST 

CONVENTIONAL  <JH  SUBJaCTIVE 

VERIFICATION 

latitude 

M2.0 

Longitude 

SM.4 

•m _ 

Fig.  C-2.  Twelve-hour-prcdiction  worksheet. 
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NORTH  AMERICAN  CYCLONE  PREDICTION  (24  HOUR) 


LATITUDE 


Ut(0) 

,  x  1.0000 

Long(i) 

‘IQ;  A- 

X  -0.0787 

P(l 1 ,9) 

iQlfc- 

X  -0.0493 

4P(n,7) 

-a.  . 

X  -0.0971 

4P(9,1) 

«=J _ 

X  0.1026 

4(7,7) 

539 

X  0.0*54 

4(9,7) 

5  35 

X  .0.2244 

4(13,7) 

ZH1. 

X  0.0863 

4(13,5) 

■SJLS- 

X  0.1188 

4(15,5) 

551. 

X  0.1476 

44(7,7) 

-  0 

X  -0.140* 

44(13,3) 

__l _ 

X  0.0883 

SUM  OF  PRODUCTS  J3.5gJ? 

LONGITUDE 

P(l3,7)  toil  X  0.0914 

P(13.1)  ion  i  -0.0800 

4P(7,5)  H  X  -0.1322 

Z(l5,t1 )  SIS  X  O.O305 

4(11,9)  53.7  X  0.1040 

4(11,5)  5W9  X  0.1070 

4(11,3)  Si 3  X  -0.2093 

4(9,1)  573.  X  0.0020 

44(11,7)  -I  X  0.1 142 

44(7,3)  -I  X  -0.1325 

H(7,5)  SAk  *  -0.0787 

SUM  OF  PRODUCTS  diU 
CONSTANT  ADDITIVE  10.1768 


CONST AMT  ADDITIVE  20.5766 
FORECAST  LATITUDE 


gam  tasssm 


LongU) 

3a.a. 

X  0.1737 

P(5,5) 

tail.- 

X  -0.1437 

P(10,5) 

lflQI- 

X  0.2071 

P(1 1 ,5) 

1CQ7- 

X  0.5390 

4P(9,5) 

x  0.4003 

4P(9,3) 

_^J _ 

X  0.1026 

4(11,1) 

57*1. 

X  -0.1217 

44(9,9) 

_ Q _ 

X  0.4669 

44(9,5) 

--7  ■ 

X  0.2011 

H( 9,7 ) 

SIH 

X  0.2566 

H(11,5) 

SH3 

X  .0.2205 

SUM  OF  PRODUCTS  542.03 


(1)  Total  ~9< _ 

(2)  1/2  (Initial  «•  fe»t.  lat.)  93.3 

(3)  Sac  (2)  1. 367  X  H)  -  -ll.SO 

(4)  Initial  Lon*,  (a)  90  3. 

(5)  FORECAST  LOW.  £(■))  +  (a)  7  7?7 


CONSTANT  ADDITIVE  434.380 


F CJT  CENTRAL  PRESSURE  _3,H, 


FORECAST 

CONVENTIONAL  OR,  SUBJECTIVE 

VERIFICATION 

Latitude 

H4.1 

Longitude 

77.7 

Contra!  Treasure 

_ 154 _ 

Fig.  C-3.  Twenty-four-hour  prediction  worksheet. 
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NORTH  AM  URIC  AN  CYCLONE  PREDICTION  (.‘Hi  HOUR) 


LATITUDE 

LONGITUDE 

ut(0) 

222- 

,  X  1.0000 

P(  1 3,7) 

iQlfc 

X  0.16*7 

Lon«U) 

90,1 

X  -0.1206 

P(t3,i) 

1017 

X  -0.1587 

P(9,7) 

iaia_ 

X  0.1168 

2(13,11) 

sal 

x  -0.0931 

P( 1 1 ,7) 

ipia 

X  -0.2193 

2(15, tl) 

Si£. 

X  0.0896 

P(  1 5,7 ) 

_1£LS_ 

X  0.0619 

2(11,9) 

S22 

X  0.1680 

^(9,7) 

'SR'S 

X  .0.2218 

2(11,3) 

Ski. 

X  -0.139* 

-05,7) 

SH5 

X  0.0666 

2(9,1) 

,  SJ2„ 

X  -0.0682 

2(13,5) 

£55 

X  0.1991 

*2(7,5) 

7.I-. 

X  -0.2001 

2(15,5) 

S-S  8 

X  0.0289 

H(7,5) 

S36 

X  -0.0820 

*2(7,7) 

_ 2 _ 

X  -0.1886 

SUM  OP  PRODUCTS  -3fe.9ff 

*H(11,5) 

—L2— 

X  0.1211 

CONSTANT  ADDITIVE  25.0906 


SUM  OP  PRODUCTS  S  S.  IH 

CONSTANT  ADDITIVE  -9.6302 

FORECAST  LATITUDE  H  S.S 

(1) 

Total  -  -H.89 

(2) 

1/2  (initial  «•  feat  lat. )  W4.Q 

CENTRAL  PRESSURE 

(3) 

Soc  (2)  1,390  X  (1)  -  ~/iS3 

LorurU)  90.  1  X  0.1900 

P( 5.7)  lOlg  X  -0.U26 

P(11.7)  .10  16  X  0.1786 

(*) 

Initial  Iona,  (a)  90.  X 

(5) 

PQtECAST  LONG,  fij)  ♦  (l)  7l3.1 

P<  10,5 )  10Q1  *  0.2955 

4P(11,7)  -a.  X  0.2990 

“’ll. 3)  Q  X  -0.1200 

-(9.9)  538  X  0.1955 

2(i3.3)  Si 7  X  .0.1127 

*(H,l)  ST-l  X  0.0)15 

*2(9,9)  Q  X  0.1937 

*2(9,5)  -1  X  0.1123 

*H(t  1 ,5)  1.7  I  -0.2962 

SUN  OP  PRODUCTS  H71.W7 
CONSTANT  ADDITIVE  523.605 
FC5T  CENTRAL  PRESSURE  995 


PORECAST 

CONVENTIONAL  OR  SUBJECTIVE 

VERIFICATION 

Latltuda 

MS.S 

Longitude 

73.7 

Cantral  Praaaura 

_ 995 

Fig.  C-4.  Thirty-six-hour-prediction  worksheet. 
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Construction  of  a  grid  overlay  similar  to  that  ill  Fig.  2-2  may  be  facilitated  by 
recalling  that  the  grid  has  the  same  spacing  between  points  as  the  JNWP  grid,  i.e., 

1  in.  on  a  polar-stereographic  projection,  with  standard  parallel  at  60*N  for  a  map  scale 
of  1:15,000,000. 

C.2  Use  ot  Worksheets 

C.2.1  Data-tabulation  Worksheet 

Using  the  grid  overlay  to  locate  the  appropriate  points,  enter  the  sea-level  pressures 
and  500- mb  heights  from  the  current  charts  (in  the  correct  units  and  format)  under  t^. 

[To  position  the  grid  overlay,  place  the  point  (k,  i)  =  (10,5)  over  the  t^-position  of  the  cyclone 
and  orient  so  that  the  line  k  =  10  coincides  with  the  meridian  passing  through  the 
^-position,  j 

Determine  the  t  entries  for  the  required  points  from  the  12-'hr-old  charts. 

[Maintain  the  gridpoint^lk./)  =  (10,5)  over  the  t^-position  of  the  cyclone;  the  values 
recorded  for  t  are  to  be  used  to  obtain  time  changes  of  heights  and  pressures  at  the 
several  points.  [ 

Determine  the  necessary  A's. 

Determine  the  1000-500-mb  thicknesses  at  the  required  points  either  by  using  the 
formulae  given  at  the  bottom  of  the  page  or  by  reading  the  values  directly  from  the 
analyzed  1000-500-mb  thickness  chart  (using  the  grid  as  an  overlay)  or  from  a  nomogram. 

C.2.2  Prediction  Worksheets 

Enter  on  each  of  the  three  prediction  worksheets  the  appropriate  parameters  from 
the  data-tabulation  worksheet. 

To  determine  the  forecast  latitudes: 

(a)  Compute  the  products  of  the  individual  parameters  and  their  respective 
coefficients. 

(b)  Sum  these  products. 

(c)  Add  the  "constant  additive"  to  this  sum;  the  total  thus  derived  is,  in  each  case, 
the  forecast  latitude. 

To  determine  the  forecast  longitudes; 

(a)  Compute  the  products  of  the  individual  parameters  and  their  respective 
coefficients. 

(b)  Sum  these  products . 

(c)  Add  the  "constant  additive"  to  this  sum;  enter  the  total  thus  derived  on  line  1 
on  the  right-hand  side  of  the  page.  [This  is  the  forecast  westward  (+)  or  eastward  (-) 
movement  in  degrees  of  latitude. ) 
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(d)  To  convert  the  forecast  movement -to  degrees  of  longitude,  take  otte*Oialf  the 
sum  of  the  initial  latitude  (0)  and  the  forecast  latitude  (determined;  aboviey.1;  This,  gives 
a  mean  latitude. 

(e)  Using  a  secant  table,  determine  the  secant  ("sec")  of  this  mean  latitude. 

Enter  the  value  on  line  3  on  the  right-hand  side  of  the  page  in  the  blank  following  thewbrds 
"sec  (2)." 

(f)  Multiply  this  secant  value  by  the  value  entered  on  line  1,  i.e.,  by  the  total  cited 
in  (c)  above.  Enter  this  product  at  the  far  right  of  line  3. 

(g)  Enter  on  line  4  the  initial  longitude  of  the  cyclone  at  time  t^.  Note  that  east 
longitude  is  entered  with  a  minus  sign,  west  longitude  with  a  plus. 

(h)  Algebraically  add  together  the  values  on  lines  3  and  4,  and  enter  this  total  on 
line  5;  this  is  the  forecast  longitude.  (Thus,  the  value  on  line  3  is  a  forecast  longitude 
change;  when  added  to  the  original  longitude,  the  forecast  longitude  is  obtained.) 

To  forecast  the  central  pressure  of  the  cyclone,  use  the  same  procedure,  but 
with  values  of  pressure. 

Unusual  displacements  or  changes  in  pressure  should  be  recomputed. 

Enter  the  forecast  latitude,  longitude,  and  central  pressure,  the  conventional  or  sub¬ 
jective  forecasts,  and  the  verification  in  the  box  at  the  bottom  of  the  page. 
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